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Abstract 

We discuss two inter-related topics: a multi-component s- and t-channel model of 'soft' 
high-energy pp interactions and the properties of the exclusive Higgs signal at the LHC. 

1 Introduction 

We begin by drawing attention to the exciting possibility of studying the Higgs sector via the 
exclusive process pp — >■ p + H + p at the LHC, where the + signs denote the presence of large 
rapidity gaps. The prediction of the event rate of such a process depends on an interesting 
mixture of 'soft' and 'hard' physics. We explain why the former requires the development of 
a multi-component s- and t-channel model of high-energy 'soft' processes, in which absorptive 
effects play a key role. We describe how the model may be used to estimate the survival 
probability of the large rapidity gaps to eikonal and enhanced soft rescattering. We comment 
on other models used to calculate the survival factors. 

We note that CDF experiments at the Tevatron have already measured the rate of similar 
exclusive processes, namely pp — > p + A + p where A = 77 or dijet or Xc PP- These processes 
are driven by the same theoretical mechanism used to estimate the exclusive Higgs signal. The 
agreement of the CDF experimental rates with the model predictions leads to optimism of the 
use of very forward proton taggers to explore the Higgs sector at the LHC [2\ . 

The discussion here is brief, with a minimum of references. More details, and references, 
can be found in two recent reviews covering the same material [31 H] . 

2 Advantages of the exclusive Higgs signal with H — >• bb 

The exclusive process pp — >■ p + H + p for the production of a Higgs at the LHC with mass 
Mn ~ 140 GeV, where the dominant decay mode is H — > bb, has the following advantages: 

• The mass of the Higgs boson (and in some cases the width) can be measured with high 
accuracy (with mass resolution o~(M) ~ 1 GeV) by measuring the missing mass to the 
forward outgoing protons, provided that they can be accurately tagged some 400 m from 
the interaction point. 

1 Based on a talk by Alan Martin at the UCL Workshop on "Standard Model discoveries with early LHC 
data" , 30 March - 1 April, 2009 
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• It offers a unique chance to study H — > bb, since the leading order bb QCD background 
is suppressed by the P-even J z = selection rule, where the z axis is along the direction 
of the proton beam. Indeed, at LO, this background vanishes in the limit of massless b 
quarks and forward outgoing protons. Moreover, a measurement of the mass of the decay 
products must match the 'missing mass' measurement. For a SM Higgs the signal-to- 
background ratio S/B ~ 0(1) 

• The quantum numbers of the central object (in particular, the C- and P-parities) can be 
analysed by studying the azimuthal angle distribution of the tagged protons. Due to the 
selection rules, the production of ++ states is strongly favoured. 

• There is a clean environment for the exclusive process — this is even possible with over- 
lapping interactions (pile-up) using fast timing detectors with very good resolution: 10 
ps or better. 

• For SUSY Higgs there are regions of SUSY parameter space were the signal is enhanced 
by a factor of 10 or more, while the background remains unaltered. Moreover, there are 
domains of parameter space where Higgs boson production via the conventional inclusive 
processes is suppressed whereas the exclusive signal is enhanced, and even such, that both 
the h and H ++ bosons may be detected. 



3 Is the exclusive Higgs cross section large enough? 

What is the price that we pay for the large rapidity gaps? How do we calculate the cross section 
for the exclusive process pp — » p + H + pi The calculation of the exclusive production of a 
heavy system is an interesting mixture of soft and hard QCD effects. The basic mechanism is 
shown in Fig. [TJ The t-integrated cross section is of the form 



S 2 

a — — 
P 2 



N J ^f g (x u x[,Ql» 2 )f g (x 2 ,a/ 2 ,QliJ, 2 ) 



(1) 



where P/2 is the t-slope of the proton- Pomeron vertex, and N is given in terms of the H — >■ gg 
decay width. The probability amplitudes, f g , to find the appropriate pairs of i-channel gluons 
(xi,x[) and (x2,x' 2 ), are given by the skewed unintegrated gluon densities at a hard scale 
fi ~ Since (x' ~ Qt/y/s) "C (x ~ M H /\fs) <C 1, it is possible to express f g (x, x', Q 2 , /i 2 ), 

to single log accuracy, in terms of the conventional integrated density g(x), together with a 
known Sudakov suppression factor T, which ensures that the active gluons do not radiate in 
the evolution from Q t up to the hard scale /i ~ Mh/2, and so preserves the rapidity gaps. The 
factor T ensures that the integral is infrared stable, and may be calculated by perturbative 
QCD. 

If we were to neglect the rapidity gaps survival factor, S 2 in ([I]), then QCD predicts that 
the exclusive cross for producing a SM Higgs of mass 120 GeV would be more than 100 fb at 
an LHC energy of y/s =14 TeV. 

The factor S 2 in (TjQ) is the probability that the secondaries, which are produced by soft 
rescattering, do not populate the rapidity gaps. As written, the cross section assumes soft-hard 
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Figure 1: The mechanism for the exclusive process pp — >■ p + H + p, with the eikonal and 
enhanced survival factors shown symbolically. The thick lines on the Pomeron ladders, either 
side of the subprocess (gg — >■ H), indicate the rapidity interval Ay where enhanced absorption 
is not permitted, see Section [7731 

factorization. In other words, the survival factor, denoted by in Fig. [H and calculated from 
an eikonal model of soft interactions, does not depend on the structure of the perturbative QCD 
amplitude embraced by the modulus signs in (0Q). Actually the situation is more complicated. 
There is the possibility of enhanced rescattering which involves intermediate partons, and which 
breaks soft-hard factorization. To calculate the corresponding survival factors, and S en h, 
we need, first, a model for soft high-energy pp interactions. We come back to the estimation of 
cr(pp p + H + p) in Section 17.31 

4 Requirements of a model of soft interactions 

Besides the need for calculating the rapidity gap survival factors, it is valuable to revisit the 'soft' 
domain at this time because of the intrinsic interest in obtaining a reliable self-consistent model 
of high-energy soft interactions which may soon be illuminated by data from the LHC Moreover, 
we need a reliable model so as to be able to predict the gross features of soft interactions; in 
particular to understand the structure of the underlying events at the LHC. 

What are the requirements of such a high-energy model? It should be self-consistent theo- 
retically - it should satisfy unitarity; absorptive corrections are large and imply the importance 
of multi-Pomeron contributions. The model should describe all the available soft data in the 
CERN-ISR to Tevatron energy range. Finally, the model should include Pomeron components 
of different size so that we can allow for the effects of soft-hard factorization breaking. 

The total and elastic proton-proton cross sections are usually described in terms of an 
eikonal model, which automatically satisfies s-channel elastic unitarity. The unitarity relation 
is diagonal in impact parameter b, and so these reactions can be described in terms of the 
opacity fl(s, b) > 

da tot /d 2 b = 2(1 - e- Q / 2 ), da el /d 2 b = (1 - e^ 2 ) 2 , (2) 

see Fig. [5^a). The Good- Walker formalism [5] is used to account for the possibility of excitation 
of the initial proton, that is for two-particle intermediate states with the proton replaced by N* 
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(a) Elastic amplitude 
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(b) Incfuslon of iow-mass dissociation 
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(c) inclusion of high-mass dissociation 
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Figure 2: (a) The single-channel eikonal description of elastic scattering; (b) the multichannel 
eikonal formula which allows for low-mass proton dissociations in terms of diffractive eigenstates 
\4>i), \4>k)', and (c) the inclusion of the multi-Pomeron-Pomeron diagrams which allow for high- 
mass dissociation. 



resonances, Fig. E(b). Diffractive eigenstates \<f>i) are introduced which only undergo 'elastic' 
scattering. That is, we go from a single elastic channel to a multi-channel eikonal, VL^. Already 
at Tevatron energies the absorptive correction to the elastic amplitude, due to elastic eikonal 
rescattering, gives about a 20% reduction of simple one-Pomeron exchange. After accounting 
for low-mass proton excitations, the correction becomes twice larger (that is, up to a 40% 
reduction). 

At first sight, by enlarging the number of eigenstates \<pi) it seems we may even allow for 
high-mass proton dissociation. However, here, we face the problem of double counting when 
the partons originating from dissociation of the beam and 'target' initial protons overlap in 
rapidities. For this reason, high-mass (M) dissociation is usually described by "enhanced" 
multi-Pomeron diagrams. The first, and simplest, such contribution to single proton disso- 
ciation do"sD / dM 2 , is the triple-Pomeron graph, see Fig. [2(c). The absorptive effects in the 
triple- Regge domain are expected to be quite large ( < 80%), since there is an extra factor of 
2 from the AGK cutting rules [BJ. Recent triple- Regge analyses [7J, which include screening 
effects, of the available data find that the bare triple-Pomeron coupling is indeed much larger 
than the (effective) value found in the original (unscreened) analyses. This can be anticipated 
by simply noting that since the original triple-Regge analyses did not include absorptive cor- 
rections, the resulting triple-Regge couplings must be regarded, not as bare vertices, but as 
effective couplings embodying the absorptive effects. That is, 



effective 

9zp — 



bare 



^ </3P 



(3) 



where S 2 is the survival probability of the rapidity gap. Due to the large bare triple-Pomeron 
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Figure 3: (a) The ladder structure of the triple- Pomeron amplitude between diffractive eigen- 
states \4>i)i \4>k) of the proton; the rapidity y spans an interval to Y = Ins. (b) A multi- 
Pomeron diagram. 



coupling (g 3 p = \g N with A ~ 0.25, where g N is the Pomeron-proton coupling), we need a model 
of soft high-energy processes which includes multi-Pomeron interactions, see, for example, the 
final diagrams in Fig. [2(c). 

5 Multi-component s- and ^-channel model 

Here we follow a partonic approach to obtain a model high-energy soft interactions [Hj. While 
the eikonal formalism describes the rescattering of the incoming fast particles, the enhanced 
multi-Pomeron diagrams represent the rescattering of the intermediate partons in the ladder 
(Feynman diagram) which describes the Pomeron-exchange amplitude. We refer to Fig. [3j 
The multi-Pomeron effects are included by the following equation describing the evolution in 
rapidity y of the opacity starting from the 'target' diffractive eigenstate \<pk)'- 

dVkM = e _ AfW , 6 )/2 e -xn k (y,b)/2 ( A + a ,d?\ Q , h) (4) 
ay \ d 2 b J 

where y' = Ins—y. Let us explain the meanings of the three factors on the right-hand-side of (j4|). 
If only the last factor, (...)Ofe, is present then the evolution generates the ladder-type structure 
of the bare Pomeron exchange amplitude, where the Pomeron trajectory otp = 1 + A + a't. 
The inclusion of the preceding factor allows for rescatterings of an intermediate parton c with 
the "target" proton k; Fig. [3^a) shows the simplest (single) rescattering which generates the 
triple-Pomeron diagram. Finally, the first factor allows for rescatterings with the beam i. In 
this way the absorptive effects generated by all multi-Pomeron diagrams are included, like the 
one shown in Fig. E^b). There is an analogous equation for the evolution in rapidity y' of 
Qi(y',b) starting from the 'beam' diffractive eigenstate \<pi). The two equations may be solved 
iteratively. 

As we are dealing with elastic amplitudes we use e~ xn ^ 2 and not e~ xn . The coefficient A in 
the exponents arises since parton c will have a different absorption cross section from that of 
the diffractive eigenstates. Naively, we may assume that the states i, k contains a number 1/A 
of partons. The factors e~ xn ^ 2 generate multi-Pomeron vertices of the form 

g « = n m X n+m - 2 g N /2 for n + m > 3 , (5) 
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where a factor 1/n!, which comes from the expansion of the exponent, accounts for the identity 
of the Pomerons. The factors n(m) allow for the n(m) possibilities to select the Pomeron Q^Q^) 
which enters the evolution OU) from the n(m) identical Pomerons. In principle, the vertices 
are unknown. However, the above ansatz is physically motivated and is certainly better than 
to assume only a triple-Pomeron coupling, that is, to assume that = for n + m > 3. 

Even though A ~ 0.2 — 0.25, the role of factors e~ XQ ^ 2 is not negligible, since the suppression 
effect is accumulated throughout the evolution. For instance, if A <C 1 the full absorptive 
correction is given by the product XQY/2, where the small value of A is compensated by the 
large rapidity interval Y. 

So far, we have allowed multi-components in the s-channel via a multichannel eikonal. 
However, a novel feature of the model of Ref. [5] is that four different t-channel states are 
included. One for the secondary Reggeon (R) trajectory, and three Pomeron states (Pi, P2, P3) 
to mimic the BFKL diffusion in the logarithm of parton transverse momentum, ln(fc t ). Recall 
that the BFKL Pomeron is not a pole in the complex j-plane, but a branch cut. Here the 
cut is approximated by three t-channel states of a different size. The typical values of k t are 
k tl ~ 0.5 GeV, k t 2 ~ 1.5 GeV and k t3 ~ 5 GeV for the large-, intermediate- and small-size 
components of the Pomeron, respectively. Thus (jlj) is rewritten as a four- dimensional matrix 
equation for in t-channel space (a = Pi, P2, P3, R), as well as being a three-channel eikonal 
in diffractive eigenstate \(pk) space. The transition terms, added to the equations, which couple 
the different t-channel components, are fixed by the properties of the BFKL equation. So, in 
principle, we have the possibility to explore the matching of the soft Pomeron (approximated by 
the large-size component Pi) to the QCD Pomeron (approximated by the small-size component 
P3). The key parameters which drive the evolution in rapidity are the intercepts 1 + A a and 
the slopes a' a of the t-channel exchanges. 

The model is tuned to describe all the available soft data in the CERN-ISR to Tevatron 
energy range. In principle, it may be used to predict all features of soft interactions at the 
LHC. All components of the Pomeron are taken to have a bare intercept A = ctp(0) — 1 = 0.3, 
consistent with resummed NLL BFKL. However, the large-size Pomeron component is heavily 
screened by the effect of 'enhanced' multi- Pomeron diagrams, that is, by high-mass dissociation, 
which results in A e g ~ 0.08 and a' eS ~ 0.25. This leads, among other things, to the saturation 
of the particle multiplicity at low pt, and to a slow growth of the total cross section. Indeed, 
the model predicts a relatively low total cross section at the LHC - a tot (LHC) ~ 90 mb. On 
the other hand, the small-size component of the Pomeron is weakly screened, leading to an 
anticipated growth of the particle multiplicity at large p t (~ 5 GeV) at the LHC. Thus the 
model has the possibility to embody a smooth matching of the perturbative QCD Pomeron to 
the 'soft' Pomeron. 

6 Long-range rapidity correlations 

We emphasize that each multi-Pomeron exchange diagram describes simultaneously a few dif- 
ferent processes. The famous AGK cutting rules [B] gives the relation between the different 
subprocesses originating from the same diagram. 
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Note that the eikonal model predicts a long-range correlation between the secondaries pro- 
duced in different rapidity intervals. Indeed, we have possibility to cut any number of Pomerons. 
Cutting n Pomerons we get an event with multiplicity n times larger than that generated by 
one Pomeron. The probability to observe a particle from a diagram where n Pomerons are cut 
is n times larger than that from the diagram with only one cut Pomeron. The observation of a 
particle at rapidity y a , say, has the effect of enlarging the relative contribution of diagrams with 
a larger number of cut Pomerons. For this reason the probability to observe another particle 
at quite a different rapidity y b becomes larger as well. This can be observed experimentally via 
the ratio of inclusive cross sections 

a mc id 2 a/dy a dy b _ d 2 N/ dy a dy b _ 

2 : (da/dy a )(da/dy b ) ~ (dN/dy a )(dN/dy b ) ' 1 j 

where dN/dy = {1/ a mc \)da / dy is the particle density. 

Without multi- Pomeron effects the value of R2 exceeds zero only when the two particles are 
close to each other, that is when the separation \y a — y b \ ~ 1 is not large. Such short-range 
correlations arise from resonance or jet production. However, multi-Pomeron exchange leads 
to a long-range correlation, R2 > 0, even for a large rapidity difference between the particles, 

\Va ~ Vb\ ~ Y. 



7 Rapidity gap survival 

Now that we have a model of high-energy soft interactions, we can estimate the rapidity gap 
survival factors S 2 ik and S 2 nh of the process pp —> p + H + p shown in Fig. [TJ We start with 



7.1 Eikonal rescattering 



The gap survival factor caused by eikonal rescattering of the diffractive eigenstates [5] , for a 
fixed impact parameter b, is 



$>,| 2 \a k \ 2 M lk {h) exp(-^° t ( S ,b)/2) 



^2\a,i\ 2 \a k \ 2 M ik (b) 

i.k 



(7) 



where ^^(s, b) is the total opacity of the ik interaction, and the a^'s occur in the decomposition 
of the proton wave function in terms of diffractive eigenstates \p) = J2i a «l0j)- The total opacity 
has the form Q^.(y)Qf (y') integrated over the impact parameters bi, b2 (keeping a fixed impact 
parameter separation b = bi — b 2 between the incoming protons) and summed over the different 
Pomeron components a. Recall y' = Y — y = Ins — y, see Fig. [3j The exact shape of the matrix 
element M.i k for the hard subprocess gg — > H in b space and the relative couplings to the 
various diffractive eigenstates i, k should be addressed further. 
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Figure 4: A "first look" at the impact parameter dependence of the signal for 120 GeV Higgs 
production at the LHC after including an eikonal rescattering correction. 

One possibility is to say that the b dependence of Ai should be, more or less, the same 
as that observed for diffractive J/ip electroproduction (7 + p — > J/ip+p), and the coupling 
to the component of the proton should be proportional to the same factor 7^ as in a soft 
interaction. This leads to 

M ik oc 7 i7fe exp(-6 2 /4fi) (8) 

with t-slope B ~ 4 GeV -2 . The resulting "first look" predictions obtained using the 'soft' 
model of [5], for the exclusive production of a scalar 120 GeV Higgs at the LHC, are shown 
in Fig. HI After we integrate over b, we find that the survival probability of the rapidity gaps 
in pp — >• p + H + p to eikonal rescattering is (S' 2 ik )=0.017, with the Higgs signal concentrated 
around impact parameter b = 0.8 fm. Expressing the survival factors in this manner is too 
simplistic and even sometimes misleading, for the reasons we shall explain below; nevertheless 
these numbers are frequently used as a reference point. 

7.2 Enhanced rescattering 

As indicated in Fig. [TJ besides eikonal screening, S^, caused by soft interactions between the 
protons, we must also consider so-called 'enhanced' rescattering, S en h, which involves interme- 
diate partons. Since we have to multiply the probabilities of absorption on each individual 
intermediate parton, the final effect is enhanced by the large multiplicity of intermediate par- 
tons. Unlike 5* 2 ik (b), the enhanced survival factor Sf^b) cannot be considered simply as an 
overall multiplicative factor. The probability of interaction with a given intermediate parton 
depends on its position in configuration space; that is, on its impact parameter b and its mo- 
mentum k t . This means that S en h simultaneously changes the distribution of the active partons 
which finally participate in the hard subprocess. It breaks the soft-hard factorization of ([T]). 

Do we anticipate that S^h will be important? Working at LO (of the collinear approxima- 
tion) we would expect that effect may be neglected. Due to strong /c r ordering the transverse 




9 



momenta of all the intermediate partons are very large (i.e. the transverse size of the Pomeron 
is very small) and therefore the absorptive effects are negligible. Nevertheless, this may be not 
true at a very low x, say x ~ 10~ 6 , where the parton densities become close to saturation and 
the small value of the absorptive cross section is compensated by the large value of the parton 
density. Indeed, the contribution of the first enhanced diagram, which describes the absorption 
of an intermediate parton, was estimated in the framework of the perturbative QCD in Ref. [9]. 
It turns out that it could be quite large. On the other hand, such an effect does not reveal 
itself experimentally. The absorptive corrections due to enhanced screening must increase with 
energy. This is not observed in the present data (see [10] for a more detailed discussion). One 
reason is that the gap survival factor S 2 ik already absorbs almost the whole contribution from 
the centre of the disk. The parton essentially only survives eikonal rescattering on the periph- 
ery; that is, at large impact parameters b. On the other hand, on the periphery, the parton 
density is rather small and the probability of enhanced absorption is not large. 



7.3 Gap survival for exclusive Higgs production 

Now, the model of Ref. [8], with its mult i- component Pomeron, allows us to calculate the 
survival probability of the rapidity gaps, to both eikonal and enhanced rescattering. Recall that 
the evolution equations in rapidity (like fll])) have a matrix form in aa! space, where a = 1, 2, 3 
correspond to the large-, intermediate- and small-size components of the Pomeron. We start the 
evolution from the large component Pi, and since the evolution equations allow for a transition 
from one component to another (corresponding to BFKL diffusion in lnk t space), we determine 
how the enhanced absorption will affect the high-/^ distribution in the small-size component 
P 3 , which contains the active gluon involved in forming the Higgs. Moreover, at each step of 
the evolution the equations include absorptive factors of the form e~ x ( Q k+ n i)/ 2 . By solving the 
equations with and without these suppression factors, we could quantify the effect of enhanced 
absorption. However, there are some subtle issues here. First, since we no longer have soft-hard 
factorization, we must first specify exactly what is included in the bare hard amplitude. 

Another relevant observation is that the phenomenologically determined generalised gluon 
distributions are usually taken at p t = 0, and then the observed "total" cross section is cal- 
culated by integrating over p t of the recoil protons assuming an exponential behaviour e -Sp *; 
that is 

f da , 2 ,2 1 da 

a = / ,9,9 d Plt d P2 



, d P lMt uplMt b 2 dplMt 

where 

'dfie-** = l/B = (p 2 ). (10) 



(9) 

Plt=P2t=0 



However, the total soft absorptive effect changes the p t distribution in comparison to that for 
the bare cross section determined from perturbative QCD. Moreover, the correct pt dependence 
of the matrix element A4 of the hard gg — » H subprocess does not have an exponential form. 
Thus the additional factor introduced by the soft interactions is not just the gap survival S 2 , 
but rather S 2 {p 2 ) 2 , where the square arises since we have to integrate over the p t distributions 
of two outgoing protons. Indeed in all the previous calculations the soft prefactor had the form 
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S 2 /B 2 . Note that, using the model of Ref. [8], we no longer have to assume an exponential b 
behaviour of the matrix element. Now the b dependence of Ad(b) is driven by the opacities, 
and so is known. Thus we present the final result in the form S 2 {p 2 ) 2 . That is, we replace S 2 / B 2 
in ([I]) by S 2 (p 2 ) 2 . So if we wish to compare the improved treatment with previous predictions 
obtained assuming B = 4 GeV -2 we need to introduce the "renormalisation" factor ((p 2 )B) 2 . 
The resulting (effective) value is denoted by S 2 S . 

Before we do this, there is yet another effect that we must include. We have to allow for 
a threshold in rapidity [10]. The evolution equation for Q, and the analogous one for 
flf, are written in the leading ln(l/x) approximation, without any rapidity threshold. The 
emitted parton, and correspondingly the next rescattering, is allowed to occur just after the 
previous step. On the other hand, it is known that a pure kinematical t min effect suppresses 
the probability to produce two partons close to each other. Moreover, this t m i n effect becomes 
especially important near the production vertex of the heavy object. It is, therefore, reasonable 
to introduce some threshold rapidity gap, Ay, and to compute S 2 nh only allowing for absorption 
outside this threshold interval, as indicated in Fig. [TJ For exclusive Higgs boson production at 
the LHC, the model gives S 2 S = 0.004, 0.009 and 0.015 for Ay = 0, 1.5 and 2.3 respectively 
For Ay = 2.3 all the NLL BFKL corrections may be reproduced by the threshold effect. 

Furthermore, Ref. [11] presents arguments that 

(S 2 S ) = 0.015 t°;°05 (11) 

should be regarded as a conservative (lower) limit for the gap survival probability in the ex- 
clusive production of a SM Higgs boson of mass 120 GeV at the LHC energy of y/s = 14 TeV. 
Recall that this effective value should be compared with S 2 obtained using the exponential 
slope B = 4 GeV -2 . The resulting value for the cross section is, conservatively, 

cr{pp^p + H + p) ~ 2-3 fb, (12) 

with an uncertainty^] of a factor of 3 up or down. 

7.4 Comments on other estimates of S 2 

A very small value S 2 nh = 0.063 is claimed in [12], which would translate into an extremely 
small value of S 2 S = 0.0235 x 0.063 = 0.0015. There are many reasons why this estimate is 
invalid. In this model the two-particle irreducible amplitude depends on the impact parameter b 
only through the form factors of the incoming protons. The enhanced absorptive effects (which 
result from the sum of the enhanced diagrams) are the same at any value of b. Therefore, the 
enhanced screening effect does not depend on the initial parton density at a particular impact 
parameter b, and does not account for the fact that at the periphery of the proton, from where 
the main contribution comes (after the S^ik suppression), the parton density is much smaller 
than that in the centre. For this reason the claimed value of S 2 nh is much too small. Besides 

2 Bcsides the uncertainty arising from that on S^ s , the other main contribution to the error comes from that 
on the unintegrated gluon distributions, f g , which enter to the fourth power. 
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this lack of kf ■<->• b correlation, the model has no diagrams with odd powers of g^p. For example, 
the lowest triple- Pomeron diagram is missing. That is, the approach does not contain the first, 
and most important at the lower energies, absorptive correction. Next, recall that in a theory 
which contains the triple-Pomeron coupling only, without the four- Pomeron term (and/or more 
complicated multi- Pomeron vertices), the total cross section decreases at high energies. On the 
other hand, the approximation used in [T2] leads to saturation (that is, to a constant cross 
section) at very high energies. In other words, the approach is not valid at high energies^. This 
means that such an approximation can only be justified in a limited energy interval; at very 
high energies it is inconsistent with asymptotics, while at relatively low energies the first term, 
proportional to the first power of the triple-Pomeron coupling g%p, is missing. Finally, the 
predictions of the model of [12] have not been compared to the CDF exclusive data of Section 

El 

The values of x ~ 10~ 6 relevant to the evaluation of S 2 nh at LHC energies are quite small. 
At leading order (LO) the gluon density increases with 1/x. So, at first sight, it appears that 
this may lead to the Black Disk Regime where the enhanced absorptive corrections will cause 
the cross section for exclusive production to practically vanish. This is the basis of the low 
value of S 2 nh estimated in [13J. However, NLO global parton analyses show that at relatively 
low scales {k 2 = fi 2 ~ 2 — 4 GeV 2 ) the gluons are flat for x < 10 -3 — 1CT 4 (or even decrease 
when x — > 0). Recall that the contribution of enhanced diagrams from larger scales decrease as 
1/k 2 (see [9]). The anomalous dimension 7 < 1/2 is not large and so the growth of the gluon 
density, xg(x,fi 2 ) oc (/i 2 ) 7 , cannot compensate the factor 1/k 2 . Therefore the whole enhanced 
diagram contribution should be evaluated at rather low scales where the NLO gluon is flatB in 
x. Thus there is no reason to expect a strong energy dependence of S 2 nh . 

Another observation against the extremely small gap survival factor, S e nh, estimated in 
[T2"l [13] is that the analysis of [11] shows the following hierarchy of the size of the gap survival 
factor to enhanced rescattering 

S^ h c (M H > 120 GeV) > SjT mn (ir,E T >5GeV) > SjT mn (Xc), (13) 

which reflects the size of the various rapidity gaps of the different exclusive processes. The fact 
that 77 and Xc events have been observed at the Tevatron confirms that there is no danger that 
enhanced absorption will strongly reduce the exclusive SM Higgs signal at the LHC energy. 



8 Exclusive processes at the Tevatron and the LHC 

In preparation for the studies of the exclusive process, pp —> p+H+p, at the LHC, it is important 
to know how reliable is the 'Durham' prediction of the cross section f[T2"j) . particularly as some 
other estimates are much lower, see Section 17.41 

Exclusive diffractive processes of the type pp — > p + A + p have already been observed by 
CDF at the Tevatron, where A = 77 [T3] or dijet [15] or Xc |16| . As the sketches in Fig. [5] show, 

3 The model of Ref. |12j should contain a parameter which specifies the energy interval where the assumed 
approximation is valid. This has not been given. 

4 Note that the flat x-behaviour of NLO gluons allow the justification of the inequalities in (TlUf below. 
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Figure 5: The mechanism for the exclusive processes observed by CDF at the Tevatron. The 
survival probabilities of the rapidity gaps are not indicated in the sketches. 

these processes are driven by the same mechanism as that for exclusive Higgs production, but 
have much larger cross sections. They therefore serve as "standard candles" . 

CDF observe three candidate events for pp — » p + 77 + p with Ej, > 5 GeV and |^ 7 | < 1 
[14J. Two events clearly favour the 77 hypothesis and the third is likely to be of 7r°7r° origin. 
The predicted number of events for these experimental cuts is O.8I0J [IZJ, giving support to 
the 'Durham' approach used for the calculation of the cross sections for exclusive processes. 

Especially important are the recent CDF data [15] on exclusive production of a pair of high 
E T jets, pp — > p + jj + p. Such measurements could provide an effective gg pp 'luminosity 
monitor' for the kinematical region appropriate for Higgs production. The corresponding cross 

a(E T > £™ in ; R j:j > 0.8) pb 

10 3 
10 2 

10 

1 

-1 

10 10 20 30 jet E™ n 

Figure 6: The cross section for 'exclusive' dijet production at the Tevatron as a function E™ m 
as measured by CDF [15]. The data integrated over the domain Rjj = M di j et /M PP > 0.8 and 
E T > E™ m . A jet cone of R < 0.7 is used. The curves are the pure exclusive cross section 
predicted by the 'Durham' model using the CDF event selection. The solid curve is obtained 
[3] by rescaling the parton (gluon) transverse momentum p p to the measured jet transverse 
energy E? by Ej- = 0.8py. The dashed curve assumes Ej- = 0.75px- The rescaling procedure 
effectively accounts for the hadronization and radiative effects, and for jet energy losses outside 
the selected jet cone. 
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section was evaluated to be about 10 4 times larger than that for the exclusive production of a 
SM Higgs boson. Since the exclusive dijet cross section is rather large, this process appears to be 
an ideal 'standard candle'. A comparison of the data with analytical predictions, obtained using 
the 'Durham' model, is given in Fig. [6j It shows the E™ m dependence for the dijet events with 
Rjj = Mdijet/Mpp > 0.8, where Mpp is the invariant energy of the incoming Pomeron-Pomeron 
system. The agreement with the theoretical expectations lends credence to the predictions for 
the exclusive Higgs production. 

Moreover, in the early data runs of the LHC it is possible to observe a range of diffractive 
processes which will illuminate the different components of the theoretical formalism. Some 
information is possible even without tagging the outgoing protons [18]. For example, the 
observation of the rapidity, y^, dependence of the ratio of diffractive (single gap) A production 
to inclusive A production will probe the effect of enhanced rescattering. The object A may be 
an T or a W boson or a dijet system. The ratio should avoid normalisation problems. Other 
informative examples are W (or Z) + rapidity gaps events or central 3-jet production. The 
exclusive process pp — > p + T + p is interesting. For low p t of the outgoing proton, the process 
is mediated by photon exchange and probes directly the unintegrated gluon distribution. At 
larger p t , the process is driven by odderon exchange and could be the first hint of the existence 
of the odderon. 

9 Conclusions 

We emphasized the value of installing near beam proton detectors to the ATLAS and CMS 
experiments in order to assist the study of the Higgs sector at the LHC via the exclusive 
process pp — >■ p + H + p. This is a unique chance to study the H — > bb decay, due to the 
large suppression of the QCD bb background, and to determine the spin, C and P values of 
the Higgs. We described how the prediction of the pp — >■ p + H + p cross section depends on a 
mixture of 'soft' and 'hard' physics. We introduced a model of 'soft' high-energy pp interactions, 
which possessed all the requirements to give a reliable estimate of the survival probability of 
the rapidity gaps to both eikonal and enhanced 'soft' rescattering effects. Finally, we noted 
that the rates of the exclusive processes already observed by CDF are in good agreement with 
the predictions of the Durham model. This lends valuable support to the exciting proposal to 
indeed install the proton taggers to explore the Higgs sector via exclusive production at the 
LHC. 
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Abstract 

We study the effect of soft gluon resummation on the gaps-between-jets cross-section 
at the LHC. We review the theoretical framework that enables one to sum logarithms of 
the hard scale over the veto scale to all orders in perturbation theory. We then present a 
study of the phenomenological impact of Coulomb gluon contributions and super-leading 
logarithms on the gaps between jets cross-section at the LHC. 

1 Jet vetoing: gaps between jets 

We consider dijet production with transverse momentum Q and a veto on the emission of ad- 
ditional radiation in the inter-jet rapidity region, Y, harder than Q . We shall refer generically 
to the "gaps between jets" process, although the veto scale is chosen to be large, Qq = 20 GeV, 
so that we can rely on perturbation theory. Thus a "gap" is simply a region of limited hadronic 
activity. 

Gaps between jets is a pure QCD process, hence the cross-section is large and studies 
can be performed with early LHC data. It is interesting because it allows one to investigate 
a remarkably diverse range of QCD phenomena. For instance, the limit of large rapidity 
separation corresponds to the limit of high partonic centre of mass energy and BFKL effects 
are expected to become important pQ. On the other hand one can study the limit of emptier 
gaps, becoming more sensitive to wide-angle soft gluon radiation. Furthermore, if one wants 
to investigate both of these limits simultaneously, then the non-forward BFKL equation enters 
the game [2]. In the following we discuss only wide-angle soft emissions. 

Accurate studies of these effects are important also in relation to other processes, in par- 
ticular the production of a Higgs boson in association with two jets. It is well known that this 
process can occur via gluon-gluon fusion and weak-boson fusion (WBF). QCD radiation in the 
inter-jet region is clearly different in the two cases and, in order to enhance the WBF channel, 
one can put a cut on emission between the jets [3J H]. This situation is very closely related to 
gaps between jets since the Higgs carries no colour charge, and QCD soft logarithms can be 
resummed using the same technique [5]. 

Given a hard scattering process, we can study how it is modified by the addition of soft ra- 
diation. If the observable is inclusive enough, then we have no effects because soft contributions 
cancel when real and virtual corrections are added together, as a result of the Bloch-Nordsieck 
theorem. However, if we restrict the real radiation to a corner of the phase space, as happens 
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for the gap cross-section, we encounter a miscancellation and are left with a logarithm of the 
ratio of the hard scale and veto scale, Q/Qo- The resummation of wide-angle soft radiation in 
the gaps between jets process was originally performed assuming that the real-virtual cancel- 
lation is perfect outside the gap, so that one needs only to consider virtual gluon corrections 
integrated over momenta for which real emissions are forbidden, i.e. over the "in gap" region 
of rapidity and with hr above the veto scale Q P [S] . We shall refer to these contributions 
as global logarithms. The resummed squared matrix element can be written as: 

\M\ 2 = l(m |e- frt e-« r |m ), 

e 2 f QdkT fh \ m 
? = - / -j—a s {k T ) , 1 

where V c is an averaging factor for initial state colour. The vector |mo) represents the Born 
amplitude and the operator T is the soft anomalous dimension: 

T = l -Yt 2 t + int a ■ t b + lp jet (y, \Ay\)(t 2 c + t 2 d ) , (2) 

where tj is the colour charge of parton i and the function pj et (Y, Ay) is related to the jet 
definition. The operator t 2 represents the colour exchanged in the i-channel: 

t 2 = (t a + t c ) 2 =t= + t= + 2t -t c . (3) 

The imaginary part of Eq. (j2J) is due to Coulomb gluon exchange. These contributions play an 
important role in the proof of QCD factorization and they are also responsible for super-leading 
logarithms [91 [10]. We notice that for processes with less than four coloured particles, such as 
deep-inelastic scattering or Drell-Yan processes, the imaginary part of the anomalous dimension 
does not contribute to the cross-section. For instance, if we consider three coloured particles, 
then colour conservation implies that t a + t b + t c = 0, and consequently 

Z7Tt a -t 6 = y (t 2 c -tl-t 2 ) , (4) 

which contributes as a pure phase. Coulomb gluons do play a role in dijet production, but they 
are not implemented in angular-ordered parton showers. We shall evaluate the impact of these 
contributions on the cross-section in the next section. 

It was later realised [TTJ that the above procedure is not enough to capture the full leading 
logarithmic behaviour. Real gluons emitted outside of the gap are forbidden to re-emit back into 
the gap and this gives rise to a new tower of logarithms, formally as important as the primary 
emission corrections, known now as non-global logarithms. The leading logarithmic accuracy 
is therefore achieved by considering all 2 — > n processes, i.e. n — 2 out-of-gap gluons, dressed 
with "in-gap" virtual corrections, and not only the virtual corrections to the 2—^2 scattering 
amplitudes. The colour structure quickly becomes intractable and, to date, calculations have 
been performed only in the large N c limit [Tj], [121 13] • 

A different approach was taken in [91 [10] , where the specific case of only one gluon emitted 
outside the gap, dressed to all orders with virtual gluons but keeping the full N c structure, 
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Figure 1: On the left we plot the gap fraction for Y = 3 (upper red curves) and Y = 5 (lower 
green curves) as a function of Q and on the right as a function of Y, for Q = 100 GeV (upper 
blue curves) and Q = 500 GeV (lower violet curves). The solid lines are the full resummation of 
global logarithms, while the dashed ones are obtained by omitting the m terms in the anomalous 
dimension. 



was considered. That calculation had a very surprising outcome, namely the discovery of a 
new class of "super-leading" logarithms (SLL), formally more important than the "leading" 
single logarithms. Their origin can be traced to a failure of the DGLAP "plus-prescription", 
when the out-of-gap gluon becomes collinear to one of the incoming partons. Real and virtual 
contributions do not cancel as one would expect and one is left with an extra logarithm. 
This miscancellation first appears at the fourth order relative to the Born cross-section and 
it is caused by the imaginary part of loop integrals, induced by Coulomb gluons. These SLL 
contributions have been recently resummed to all orders in |14] . The result takes the form: 

\Mf L | 2 = - 2 - [ Q ^a s (k T ) L f) (a? + n»n) j (5) 

where ^#jy) is the resummed real (virtual) contribution in the limit where the out-of-gap gluon 
becomes collinear to one of the incoming partons. The presence of SLL has been also confirmed 
by a fixed order calculation in pE]; in this approach SLL have been computed at 0(a b s ) relative 
to Born, i.e. going beyond the one out-of-gap gluon approximation. 

2 LHC phenomenology 

In this section we perform two different studies. Firstly we consider the resummation of global 
logarithms and we study the importance of Coulomb gluon contributions, comparing the re- 
summed results to the ones obtained with a parton shower approach. We then turn our at- 
tention to SLL and we evaluate their phenomenological relevance. In both studies we consider 
y/S = 14 TeV, Q = 20 GeV, jet radius R = 0.4 and we use the MSTW 2008 LO parton 
distributions [T6] . 
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Figure 2: The gap cross-section obtained using HerwigH — h (black histogram) is compared to 
the one from resummation (red curves). As before the solid line is the full result, while the 
dashed line is obtained by omitting the Coulomb gluon contributions. At the bottom we plot 
the ratio between the results obtained from the resummation and the one from HerwigH — h 



Soft logarithmic contributions are implemented in HerwigH- + via angular ordering of suc- 
cessive emissions. Such an approach cannot capture the contributions coming from the imagi- 
nary part of the loop integrals, due to Coulomb gluon exchange. We evaluate the importance 
of these contributions in Fig. [TJ On the left we plot the gap cross-section, normalised to the 
Born cross-section (i.e. the gap fraction), as a function of Q at two different values of Y and, 
on the right, as a function of Y at two different values of Q. The solid lines represent the re- 
sults of the resummation of global logarithms; the dashed lines are obtained by omitting the m 
terms in the soft anomalous dimension matrices. As a consequence, the gap fraction is reduced 
by 7% at Q = 100 GeV and Y = 3 and by as much as 50% at Q = 500 GeV and Y = 5. 
Large corrections from this source herald the breakdown of the parton shower approach. In 
Fig. [2] we compare the gap cross-section obtained after resummation to that obtained using 
HerwigH-H- [T7J [HI [IH [20] after parton showering (Q is taken to be the mean pr of the two 
leading jets). The broad agreement is encouraging and indicates that effects such as energy 
conservation, which is included in the Monte Carlo, are not too disruptive to the resummed 
calculation. Nevertheless, the histogram ought to be compared to the dotted curve rather than 
the solid one, because HerwigH — h does not include the Coulomb gluon contributions. The 
resummation approach and the parton shower differ in several aspects: some non-global loga- 
rithms are included in the Monte Carlo and the shower is performed in the large iV c limit. Of 
course the resummation would benefit from matching to the NLO calculation and this should 
be done before comparing to data. 

Finally we want to study the relevance of the SLL contributions. In order to do that we 
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Figure 3: On the left we plot the K- factor as defined in Eq. (jHJ) as a function of Q for Y = 3 
(upper red curve) and Y = 5 (lower green curve); on the right we plot it as a function of Y, for 
Q = 100 GeV (upper blue curve) and Q = 500 GeV (lower violet curve) 
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(1) u 

0.(0) ' W 

where contains the resummed global logarithms and the resummed SLL contribution 
coming from the case where one gluon is emitted outside of the gap. The results are shown in 
Fig. |3j Generally the effects of the SLL are modest, reaching as much as 15% only for jets with 
Q > 500 GeV and rapidity separations Y > 5. The contribution coming from n > 2 out-of-gap 
gluons is thought to be less important [H] . Remember that we have fixed the value of the veto 
scale Qq = 20 GeV and that the impact will be more pronounced if the veto scale is lowered. 



3 Conclusions and Outlook 

There is plenty of interesting QCD physics in "gaps-between-jets" and measurement can be 
performed with early LHC data. There are significant contributions arising from the exchange 
of Coulomb gluons, especially at large Q/Qq and/or large Y, which are not implemented in the 
parton shower Monte Carlos. However before comparing to data, there is a need to improve 
the resummed results by matching to the fixed order calculation. These observations will have 
an impact on jet vetoing in Higgs-plus-two-jet studies at the LHC. 

We have studied the super-leading logarithms that occur because gluon emissions that are 
collinear to one of the incoming hard partons are forbidden from radiating back into the veto 
region. Even if their phenomenological relevance is generally modest, they deserve further study 
because they are deeply connected to the fundamental ideas behind QCD factorization. 
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Jet-based approach for underlying-event characterization 



We discuss which characteristics of the underlying event might be useful to measure 
for improving understanding of its properties and simulating it well in Monte Carlo gener- 
ators. We use a method of selection of the underlying event which is based on jets analysis 
and does not require any phase space cuts. We analyze the predictions of several tunes of 
simulation programs for selected quantities. 

4 Introduction 

In pp collisions the hard interaction is always accompanied by soft background activity, called 
the underlying event (UE), involving the spectator partons. Underlying event has a largely 
non-perturbative nature and its separation from the hard process is to some extent a matter 
of definition. Therefore, our understanding of UE is rather incomplete and it is not clear how 
the UE should be modeled. Yet, the underlying event is very important and has a big impact 
on the results of data analysis since it adds a considerable amount of transverse momentum to 
the event. 

Improper estimation of UE may lead e.g. to modification of the inclusive jet spectrum by 
up to 50% (at moderate transverse momenta), incorrect kinematic reconstruction of mass peak 
position, degradation of the mass peak as well as lowering the efficiency of lepton isolation. 
Therefore, it is important that the models implemented in Monte Carlo generators (MC) de- 
scribe various aspects of UE as well as possible. This requirement leads to the question about 
the relevant properties of UE that could help to better tune the generators as well as further 
constrain the models. Of particular value would be the characteristics for which the current 
MCs and tunes give different predictions when extrapolated to the LHC energy. 

In the usual method of the underlying event analysis, which we call 'topological', as a 
first step, one identifies the direction of the leading jet and defines the cone around it in the 
transverse plane. Then, the UE is measured using only those regions which are transverse with 
respect to the leading jet. This method was widely used in Tevatron data analysis [U El E]- 

In these proceedings we present the study of the underlying event using an alternative 
technique based on the method for pileup and UE estimation proposed in jl]. 
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5 The method 



The method used in our analysis is jet-based, exploits the concept of jet areas [5] and can be 
carried out using the Fast Jet package [HI [7]. It proceeds as follows: 

For each event, all true particles together with a large number of extremely soft particles 
(called 'ghosts') are clustered with an infrared and collinear safe jet finding algorithm. This 
way, we obtain a list of jets with transverse momenta, p t , ranging from the scale of hard collision 
to essentially zero (for the so called 'pure ghost jets'). The susceptibility to the contamination 
from the soft radiation for each jet j is given by the jet area Aj, proportional to the number 
of ghosts in the jet. Hence, the amount of p t added by this radiation to jet j is pAj, where p 
is the level of transverse momentum per unit area characteristic for the UE. 

In order to measure p, the median of the p t j/Aj distribution is determined for the single 
event 



The main advantage of using the median is its very small sensitivity to perturbative radiation, 
which starts to be relevant only at the order a 2 s A — a A J (depending on jet definition and rapidity 
range) [I]. At the same time, there is no need for cuts in (y, </>)-plane. This is also important 
since such cuts can potentially introduce a bias (see [8]). Since the separation between UE and 
the hard part of the event is based solely on the value of ptj/Aj we refer to this method as a 
'dynamical' selection. 

One is also interested in the uncertainty of p caused by fluctuations of UE from point to point 
within the event. This quantity, denoted by a, is determined from the sorted list of {p t j/Aj} 
and given by the value for which 15.86% of jets have smaller p t j/Aj. With such definition, in the 
case of Gaussian distribution of UE, 68.27% of jets satisfy p — a / \I~Aj < p t j/Aj < p + a / \f~Aj. 




6 The results 

In this study we are interested in the UE itself. We address the question which properties of 
UE are relevant and what are the predictions for these properties from different Monte Carlo 
tunes. One important extension with respect to [H [9] is that we go differential in rapidity. 
We note that the dynamical method works well only if the number of jets used to calculate 
median is large enough. This sets some limit on the smallest possible rapidity range that can 
be used for p determination. 

We study the UE characteristics for a series of Monte Carlo generators/tunes: DW and DWT 
tunes of Pythia 6.4 |T0] by R. Field that come with the 'old' shower and SOA tune pH [T2~l [13] 
by P. Skands with the 'new' shower, Herwig 6.5 (HI [15] + Jimmy 4.3 [16] in an Atlas tune by 
A.Moraes [2] as well as the Herwig 6.5 default UE model. All models of UE, except for Herwig 
default case, are based on multiple parton interactions. 

We consider dijet events at the LHC, y/s = 14 TeV, generated with pt >m i n = 50 GeV. The 
analysis is performed with Fast Jet 2.4.1 [?]. Except where stated, we cluster particles using 
the Cambridge/ Aachen algorithm with R = 0.6 and active area with explicit ghosts. The 
results discussed in these proceedings were obtained without further rapidity-based selection. 




(7) 
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Figure 4: Left: Rapidity dependence of the average level of UE for various generators/tunes. 
Right: Average level of point-to-point fluctuations of UE normalized to (p) as a function of 
rapidity. 



We have established that for such a case one can safely go with size of the rapidity strip down 
to Ay = 1. For the analysis of the impact of further rapidity selection we refer to [8]. Here, 
we note only that the number of jets in a given range of rapidity of size of one unit (which is 
typically around 15) may not be sufficient if one requires e.g. the two hardest jets to lie also 
in this range. In such case, one can improve the quality of p determination either by removing 
the hardest jets or by doubling the size of the rapidity strip. 

In the left panel of Fig. @] we show the average value of p as a function of y, determined 
from 2 • 10 5 events for different tunes. We observe a significant dependence of the level of UE 
on rapidity. The predictions for the UE at the LHC differ also for different tunes. They are 
very close for Pythia DW and SOA since both tunes have the same yfs scaling of pt, mm- This 
scaling is faster than in Pythia DWT tune hence the latter gives larger level of UE, which is 
similar to the predictions of Herwig+Jimmy in Atlas tune. 

In the right panel of Fig. H] we show the level of fluctuations of the UE within an event 
from different tunes. It is given in terms of the average a normalized to the average p. We see 
that these fluctuations are large and that predictions from Pythia and Herwig+Jimmy differ 
significantly. 

But UE fluctuates also from event to event. To illustrate this, we plot in Fig. Othe distribu- 
tion of p determined in rapidity strip < y < 1. It is this distribution (and analogous ones for 
different rapidity ranges) that we used to compute (p) depicted in the left hand side of Fig. HJ 
The leftmost bin with the high yield contains very quiet events with p < 0.25 GeV per unit 
area. We notice sizable differences in shape for various generators/tunes. All distributions, 
however, are fairly broad, which indicates high fluctuations of UE from one event to another. 
This constitutes an argument in favor of an event-by-event analysis of the underlying event. 
In the right hand side of Fig. [5] we show the standard deviation (normalized to the average p) 
for the distribution from the left hand side of this figure and similar ones for different rapidity 
bins. We see, in particular, that the level of these fluctuations is considerably larger than the 
average level of fluctuations within an event shown in Fig. HJ 
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Figure 5: Left: Distribution of p determined in rapidity range < y < 1 from 2 • 10 5 generated 
events for various generators/tunes. Right: Standard deviation of p distributions like the one 
from left hand side, normalized to (p), as a function of rapidity. 

Correlations between values of p in different rapidity strips is another interesting character- 
istic of the underlying event. To quantify it, we use the standard definition of the correlation 
coefficient 

/ x {p{yi)p{y2)) - {p{yi))(p{y2)) , Q . 

corr(yi,y 2 ) = ==, (8) 

V<p(z/i) 2 > - (p(yi)) vW) 2 ) " MwO> 

where yi, y% are the rapidity bins of width Ay = 1 and the (...) denotes the average over 
many events (2 • 10 5 , in the case of this study). In Fig. [6] we show the coefficient from Eq. (jSJ) 
as a function of y<i for two bins of y\\ < y\ < 1 (left) and 3 < y% < 4 (right). We see that 
the predictions for the amount of correlations of the underlying event at the LHC are very 
different for Pythia and Herwig+Jimmy. They are fairly similar for various Pythia tunes is 
spite of (in some cases, like DWT or DW vs. SOA) important differences in UE models that 
come with them. We notice that the difference between Pythia and Herwig+Jimmy is the most 
pronounced discrepancy among all quantities that we have studied, which makes the correlation 
coefficient a particularly valuable quantity from the point of view of improving the fits as well 
as the model of UE. We also remark that the observed higher correlations for Herwig+Jimmy 
with respect to Pythia is qualitatively consistent with the lower level of average point to point 
fluctuations in the event discussed previously in the context of Fig. HI 

The results presented up to this point were obtained with R = 0.6, a choice motivated by 
the findings from [I] . We conclude our discussion of selected characteristics of UE that are of 
potential interest by examining the dependence of the level of UE, p, on the radius R in the 
jet definition. Let us first consider, for the purpose of illustration, a simple two-component 
model of UE in which pt/A of soft particles is given by the Gaussian distribution with mean po- 
The second component is just n h hard (with respect to the scale of the Gaussian component) 
particles. It can be shown [8] that for the case of this model 

p = Q(R — Rq)(po + const ■ o n h R). (9) 
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Figure 6: Correlations of the underlying event between different rapidity ranges. The coefficient 
from Eq. (JSJ) is plotted as function of 1/2 for two y\ bins and for a range of generators/tunes. 



The value of R has to be larger than some critical value Ro so that the number of jets that 
contain physical particle is larger than the number of ghost jets and the median formula (J7J) 
gives the value considerably larger than zero. Further increasing of R leads to linear growth of 
estimated p with the coefficient which is proportional to the number of hard jets as well as to 
the width of the Gaussian distribution a / a/ (A) (this behavior is observed also in a toy model 
discussed in [H]). The qualitative features of Eq. (jUj) turn out to be reproduced by the realistic 
simulation of dijet production at the LHC, as depicted in Fig. [7J which shows p, determined 
from the wide rapidity range \y\ < 5, for a typical, single event and for three different Monte 
Carlo models of UE. 

For better understanding, in Fig. [S] we show the corresponding histograms of pt/A. We 
see that in the case of Herwig default there are practically no hard jets except those from the 
dijet system and the whole noise is concentrated around a single soft scale. That is why the 
corresponding curve in Fig.[7Jis almost flat for R greater than some threshold value, in agreement 
with Eq. (Q for very small n^. The distribution of Pt/A for the case of Herwig+ Jimmy acquires 
some tail due to multiple parton interactions present in this model. Hence, the event contains 
larger number of jets with hard or semi-hard scale. This is consistent with a non-zero slope 
of p(R) for Herwig-I- Jimmy in Fig. [7J The slope for Pythia, in the same figure, is even larger 
which again is qualitatively consistent with a longer tail (more hard jets) in Fig. [HJ 



7 Conclusions 

We presented the analysis of the underlying event predictions for the LHC. We used the jet- 
based method of dynamical selection of UE introduced in [1] . In this approach the underlying 
event is analyzed on an event-by-event basis and one does not impose cuts in (y, 0)-plane, hence 



26 



Herwig default 
Herwig+Jimmy (Atlas) 
Pythia (DWT) 



Cambridge/Aachen 
|y|<5 



0.2 0.4 0.6 0.8 1 1.2 

R 



Figure 7: Dependence of the level of the underlying event, p, on the jet radius R for a single, 
typical event from different generators/tunes. 
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Figure 8: Histograms of pt/A corresponding to the events shown in Fig. [71 

the whole event is used. One can, however, analyze the UE locally in limited ranges of rapidity 
(or/and azimuthal angle) to study how it depends on the region of phase space. 

We discussed several quantities which we find useful as the UE characteristics: transverse 
momentum density per unit area, p, point-to-point fluctuations, a, event-to-event fluctuations, 
p distributions over many events, point-to-point correlations. We analyzed how these quantities 
depend on rapidity and jet radius for various Monte Carlo models and tunes. 

We found that the existing MC models predict for the LHC the amount of 1-3 GeV of 
transverse momentum per unit area coming from the underlying event. This value depends on 
rapidity, varies from one event to another, and changes with jet radius R used as a parameter 
of the jet finding algorithm. The amount of the UE predicted for the LHC differs also between 
generators/tunes. 

We observe large fluctuations of UE in the simulated data both within a single event (40- 
70%) and from event to event (90-130%). Here, the differences between MCs are even more 
pronounced than for the case of p. But the quantity which shows the largest discrepancy 
between Herwig+Jimmy and Pythia is the correlation coefficient. 
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We believe that all these makes the discussed characteristics very promising for further MC 
tunes as well as improving the models of the underlying event. 
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Abstract 

We review selected aspects of exclusive diffractive photoproduction at the Tevatron 
and discuss the prospects for the early LHC running. This talk [T] is based on the results 
presented in Ref. [2], with some updates due to recent experimental results. 

1 Introduction 

Exclusive diffractive Higgs boson production at the LHC, pp — > p + H+p, has attracted increas- 
ing attention in recent years as an alternative way to explore the Higgs sector [3]. However, 
the theoretical uncertainties are comparatively large relative to the inclusive production. In 
particular, the generalised (or skewed) unintegrated gluon distribution, which can be written 
in terms of the usual gluon distribution, enters to the fourth power, and is required in the 
region of small x and low scales ~ 1-2 GeV where the gluon distributions obtained from global 
fits have large uncertainties. While current Tevatron and early LHC data will prove vital in 
checking the predictions, it is also important to look for complementary processes to constrain 
the various ingredients of the calculation. One such process, where the gap survival factor is 
expected to be much closer to 1, is exclusive photoproduction, •yp — > E + p, where the photon 
is radiated from one of the two incoming hadrons; see Fig. [TJ 

fti(-Pi) •) hAP[) 



Figure 1: Exclusive photoproduction of a massive final state E = J/ip, T, Z/y* in hadron- 
hadron collisions. 
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2 Exclusive photoproduction 



Exclusive diffractive vector meson production, ^*>p — > V + p, and deeply virtual Compton 
scattering (DVCS), ^*p — > 7 + p, have been extensively studied at HERA. These processes 
provide a valuable probe of the gluon density at small x [H El E] . 

To obtain the hadron-hadron cross section for exclusive production of a massive final state 
E with rapidity y, we need to multiply the photon-hadron cross section by the flux dn/dk of 
quasi-real photons with energj]^ k ~ (M s /2) exp(y) ~ W 2 /(2y/s) [TJ: 

do" , , , , ,_, , . , dn . ^ 

—(h 1 h 2 ^h 1 + E + h 2 ) = k—a{jp^E + p), 1) 
dy dk 

together with a second term with y — > —y to account for the contribution from the interchange 
of the photon emitter and the target, neglecting interference. We also neglect absorptive correc- 
tions, and only present cross sections integrated over final-state momenta, then these effects are 
expected to be largely washed out, with a rapidity gap survival factor S 2 ~ 0.7-0.9. Alternative 
calculations including a detailed treatment of these effects can be found in Refs. [SI [10l [11] . 

The photon energy spectrum dn/dk in Eq. CD is given by a modified equivalent-photon 
(Weizsacker- Williams) approximation [7] [T2"]: 



dn a P 



dk 2irk 



2fc x 2 



1 + 



, , 11 3 3 1 , 

lnA -¥ + I-2^ + 3l5l' < 2 > 



where A = 1 + (0.71 GeV 2 )/(5^ in , Qmin — ^ 2 /7l an d 7l = ^/{2m p ) is the Lorentz factor of a 
single beamQ 

If the photon-induced contribution to exclusive vector meson production is known precisely 
enough, there is potential for odderon discovery. The odderon is the C-odd partner of the 
Pomeron, and in perturbative QCD is modelled by exchange of three gluons in a colour-singlet 
state. Cross sections were calculated using /cr-factorisation in Ref. [13J for different scenarios. 
However, no attempt was made to tune the photoproduction contribution to HERA data, 
hence the odderon-to-photon ratios are more reliable than the absolute cross sections. The 
odderon-to-photon ratios for da/dy\ y= Q at the Tevatron are 0.3-0.6 (J/ip) and 0.8-1.7 (T), 
while the LHC ratios are smaller at 0.06-0.15 (J/ip) and 0.16-0.38 (T) [13]. Moreover, odderon 
exchange leads to a different meson px distribution than the photon-induced contribution. In 
the following, we restrict attention to photoproduction. 

To calculate the exclusive 7J9 cross section in Eq. [T] we use the dipole model approach for 
exclusive diffractive processes [H] , where the amplitude factorises into the light-cone wave func- 
tions of the incoming and outgoing particle and a dipole cross section describing the interaction 
of the qq dipole with the proton. The parameterisation of the dipole cross section is given in 
terms of a DGLAP-evolved gluon density, fitted to F 2 data, with Gaussian impact parame- 
ter (6) dependence, denoted "b-Sat" model [El [TB] , which has already been shown to give a 



6 The photon-hadron and hadron-hadron centre-of-mass energies are labelled W and y/s, respectively. 
7 The results of Ref. [5] were erroneously obtained with A 2 instead of A 3 in the last term of Eq. [5J The results 
presented here have been corrected, although the numerical impact is completely negligible. 
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Figure 2: Exclusive J/^ photoproduction in hadron-hadron collisions. No absorptive correc- 
tions are included. 



good (parameter-free) description of a wide variety of HERA data on exclusive 7*p — > V + p 
(V = p,(j),J/ip,j) and inclusive Ff, F?, F L , F 2 D . Note that although the "b-Sat" model in- 
corporates saturation effects via the eikonalisation of the gluon density, these saturation effects 
are expected to be only moderate for J/ip production and negligible for T and Z° production. 
We use a "boosted Gaussian" vector meson wave function pTTj [T6]. 



3 Results for Tevatron and LHC 

3.1 Exclusive J/ip production 

In Fig. [2(a) we show the jp cross section as a function of W, where we indicate the W values 
corresponding to central production (y = 0) at the Tevatron and LHC. The "b-Sat" model 
predictions are normalised to best fit the HERA J/ip data [18], [19] by a factor 1.08. Also shown 
are the results of a direct power-law fit to the HERA data, which gives oi^ip — >■ J/tp + p) = 
(2.96 nb)(Wy GeV) ' 721 . In Fig. 0(b) we show the rapidity distributions at the Tevatron and 
LHC given by Eq. [[J The fact that y = corresponds to W values where precise HERA data 
are available means that the uncertainties in the predictions are small. CDF have recently 
measured da/dy\ = (3.92±0.62) nb [20]. The "b-Sat" model prediction of 3.4 nb multiplied 
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Figure 3: Exclusive T photoproduction in hadron-hadron collisions. No absorptive corrections 
are included. 

by S 2 ~ 0.9 [9] and an estimated odderon contribution of a factor 1.3-1.6 [13J gives a total theory 
prediction of (4.0-4.9) nb, in agreement with the Tevatron data. At the LHC, measurement of 
exclusive J / ij) production is unlikely to be possible by ATLAS or CMS due to lack of a low-pr 
trigger on leptons, but the measurement should be possible by ALICE [21] and by LHCb. 

3.2 Exclusive T production 

Only very sparse data are available on exclusive T photoproduction at HERA, and conse- 
quently, there is a much larger uncertainty in the theory predictions extrapolated to the 
Tevatron and LHC. Indeed, until very recently, only two data points with very large errors 
were published [221 123], together with a further two preliminary ZEUS data points [24J. A 
power-law fit made in Ref. [2] to these four data points, shown in Fig. Ela), gave a{^p — > 
T + p) = (0.119nb)(iy/GeV) 1,63 . The two preliminary ZEUS data points, especially the point 
at highest W, moved down slightly in the final measurement |25j, with a reduction in the 
size of the uncertainty, so that a revised power-law fit to the four published data points gives 
(j(7j9 -> T + p) = (0.968 nb)(W/ GeV) L14 , also shown in Fig. Ha). The "b-Sat" model pre- 
dictions have a very similar W dependence to this revised power-law fit, but lie more than a 
factor two below the data for the default choice of m& = 4.5 GeV and the "boosted Gaussian" 
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T wave function. Given that there are uncertainties in these two choices, which are expected 
to mainly affect the normalisation but not the W dependence, we have simply rescaled the 
"b-Sat" predictions by a factor 2.16 to best fit the final HERA T data. In Fig. El^b) we show 
the rapidity distributions at the Tevatron and LHC given by Eq. [IJ Note that the results of 
Ref. [2] for the LHC rapidity distribution, using the preliminary ZEUS T data, indicated a large 
difference between the rescaled "b-Sat" predictions and the HERA power-law fit. But with the 
new published ZEUS data, the two predictions are in much better agreement. Of course, the 
fact that the central value of the power-law fit changes so much when two of the data points 
shift within their experimental errors, see Fig. |3]^a), means that the uncertainty on the power- 
law parameterisation is large. Ideally, the errors on the two parameters in the power- law fit 
(given by the experimental errors on the 4 HERA data points), would be propagated through 
to the Tevatron and LHC rapidity distributions. An early attempt in this direction, but only 
for the error in the normalisation and not in the power, was made in Ref. [TJ. 

Candidate exclusive T events have been found by CDF and the cross section measure- 
ments are eagerly awaited. Note that the odderon contribution to exclusive T production at 
the Tevatron is predicted to be about the same or even greater than the photon-induced con- 
tribution [T3]. A feasibility study has been carried out by CMS for 100 pb" 1 of integrated 
luminosity [26J. 

The effect of the decay lepton acceptance of the Tevatron and LHC experiments has been 
studied in Ref. |27J, where the results of the power-law fit from Ref. [2] (using the preliminary 
ZEUS data) were compared with predictions obtained using an alternative parameterisation of 
the dipole cross section. Note that LHCb has the potential to measure exclusive T production 
for more forward lepton pseudorapidities than ATLAS or CMS. 

3.3 Exclusive Z/^f* production 

DVCS at HERA, 7*p — > 7 + p, is theoretically cleaner than exclusive vector meson production 
since there is no uncertainty from the wave function. The existing data are well described by the 
"b-Sat" dipole model EE] , although they are not as precise as the HERA data on exclusive 
J/ifj production. The DVCS process in ep scattering interferes with the purely electromagnetic 
Bethe-Heitler (BH) process where the real photon is instead emitted from either the incoming 
or outgoing electron. The BH process is precisely calculable in QED and is therefore subtracted 
in existing DVCS measurements at HERA. Analogous processes to DVCS at hadron-hadron 
colliders are exclusive Z° photoproduction, jp — > (Z° — > £ + £~) + p, and timelike Compton 
scattering (TCS), 7p — > (7* — > £ + £~) + p. Similarly to the DVCS case, there is interference 
of TCS with the pure QED subprocess (77 — > £ + £~) which is precisely calculable and can be 
reduced with suitable cuts [29J. 

Wave functions for an outgoing Z/'j* with timelike q 2 = M 2 > were derived in Ref. [2]. 
Differences were found with respect to the usual spacelike case, q 2 = —Q 2 < 0, such that 
the amplitude for ^p — > Z° + p is not simply the DVCS amplitude at Q 2 = M§ with a 
different coupling. In particular, we pick up a real contribution to the amplitude related to 
the contribution of an on-shell qq pair in addition to the usual imaginary part. In the dipole 
picture, direct numerical integration over the dipole size r proved to be difficult due to a 
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Figure 4: Exclusive Z° photoproduction in hadron-hadron collisions. No absorptive corrections 
are included. 

wildly oscillatory integrand. This problem was solved by taking the analytic continuation to 
complex r, then choosing an appropriate integration contour [2]. Alternatively, there are no 
such problems if working in transverse momentum space and using /^-factorisation [UJ. We 
show the results in Fig. HI where it can be seen that the cross sections at y — are enhanced 
by 5% at the Tevatron and 21% at the LHC in the (correct) timelike case compared to the 
(incorrect) spacelike case. (The odderon contribution to exclusive Z° production is expected 
to be strongly suppressed [2].) The recent calculations of Ref. [TT] found cross sections larger 
by a factor ~3, presumably due mainly to differences in the gluon density at the relevant 
x ~ Mz/yfs, while absorptive corrections were found to lower the cross section by a factor 
1.5-2. 

CDF have made a search for exclusive Z° production at the Tevatron |3D]- Eight candidate 
events were found in 2.20 (2.03) fb -1 of data in the electron (muon) channel with Mm > 40 
GeV and r\i < 4, consistent with the QED prediction for 77 — > . No candidate events were 
found in the Z° mass window, allowing an upper limit to be placed for the exclusive Z° cross 
section of a < 0.96 pb at 95% confidence- level, compared to the theory prediction of 0.3 fb [2], 
i.e. 3000 times lower than the experimental limit. The theory prediction of a — 13 fb [2] at the 
LHC looks slightly more promising. 
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da/dy\ y=0 (nb) 


o (nb) 


Event rate (s x ) 


Tevatron 
LHC 


3.4 
9.8 


28 
120 


0.33 
71 




T{1S) 


da/dy\ y=0 (pb) 


a (pb) 


Event rate (hr x ) 


Tevatron 
LHC 


10 

53 


83 
771 


1.5 

688 




Z° 


da/dy\ y=0 (fb) 


a(fb) 


Event rate (yr 1 ) 


Tevatron 
LHC 


0.077 
1.4 


0.30 
13 


0.065 
134 



Table 1: Event rates include leptonic branching ratio and assume a luminosity C = 2 x 
10 32 cm~ 2 s _1 (Tevatron) and C = 10 cm~ 2 s (LHC). No gap survival factor included. 

4 Summary 

A summary table of predictions is given in Table [TJ where the event rates (but not the cross 
sections) include the appropriate leptonic branching ratio. The T cross sections have been 
revised with respect to Ref. [2J. Note that the Tevatron and LHC design luminosities are 
assumed in all cases as an illustrative comparison of the event rates for different processes, but 
these are not realistic values for early LHC running. In particular, exclusive J / ip production is 
likely to be measured only at ALICE where the nominal luminosity (and so the event rate) is 
smaller by a factor of 2000 and at LHCb where the corresponding reduction factor is 50. 
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Abstract 

The TOTEM experiment at the LHC is dedicated to the measurement of the total 
proton-proton cross-section and to the study of the elastic scattering and diffractive dis- 
sociation processes. The main features of the TOTEM experimental apparatus and of its 
physics programme are here presented, together with some prospects for early measure- 
ments in the hrst year of the LHC. 

5 Introduction 

TOTEM p] foresees specific measurements and experimental techniques which are very differ- 
ent from the other 'general purpose' experiments at LHC. A precise 'luminosity independent' 
measurement of o v ^ OT ^ based on the Optical Theorem, will be achievable in special beam optics 
runs by simultaneously measuring: 1) the elastic scattering rate at low transfer momentum, 
possibly as small as \t\ = 1(C 3 GeV 2 , and 2) the inelastic scattering rate with the largest 
possible coverage to reduce losses to few percents. The first goal requires detectors located 
into units mounted into the vacuum chamber of the accelerator, called Roman Pots, as the 
scattered protons are emitted at angles of the order of 10 //rad, therefore without leaving the 
beam-pipe. The latter requires the measurement of all the inelastically produced particles in 
the very forward direction with respect to the pp collision point; this can be achieved by using 
tracking detector telescopes with a complete azimuthal coverage around the beam-pipe. A 
flexible trigger provided by TOTEM detectors will allow to take data under all LHC running 
scenarios. The combination of the CMS and TOTEM experiments will also allow the study 
of a wide range of diffractive processes with an unprecedented coverage in rapidity. For this 
purpose the TOTEM trigger and data acquisition (DAQ) systems are designed to be compat- 
ible with the CMS ones, in order to allow common data taking periods foreseen at a later 
stage [2J. Finally, the aim of the TOTEM experiment to obtain accurate information on the 
basic properties of proton-proton collisions should also provide a significant contribution to the 
understanding of very high energy cosmic ray physics. In the following, after a general overview 
of the experimental apparatus, the main features of the TOTEM physics programme will be 
described. 
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Figure 5: Top: The Tl and T2 telescopes embedded into the forward region of the CMS 
detector. Bottom: Roman Pots location along the LHC beam-line. All TOTEM detectors are 
located on both sides of IP5. 

6 Experimental Apparatus 

Located on both sides of the interaction point IP5 (shared with the CMS experiment), the 
TOTEM experimental apparatus comprises "Roman Pot" (RP) detectors and the Tl and T2 
inelastic telescopes (Fig. [5]). The RPs, placed on the beam-pipe of the outgoing beam in two 
stations at about 147 m and 220 m from IP5, are special movable beam-pipe insertions designed 
to detect "leading" protons with a scattering angle down to few /xrad. Tl and T2, embedded 
inside the forward region of CMS, provide charged track reconstruction for 3.1 < < 6.5 
(77 = — ln(tan^)) with a 2n coverage and with a very good efficiency. These detectors will 
provide a full inclusive trigger for all inelastic and diffractive events, minimizing losses to a few 
percent, and will be also used for the reconstruction of the event interaction vertex, so to reject 
background events [3]. The read-out of all TOTEM sub-detectors is based on the digital VFAT 
chip [3], specifically designed for TOTEM and characterized by trigger capabilities. 

The RPs host silicon detectors which are moved very close to the beam when it is in stable 
conditions. Each RP station is composed of two units in order to have a lever arm for local 
track reconstruction and trigger selections by track angle. Each unit consists of three pots, 
two vertical and one horizontal completing the acceptance for diffractively scattered protons. 
Each pot contains a stack of 10 planes of silicon strip detectors (Fig. El left). Each plane has 
512 strips (pitch of 66 /im) allowing a single hit resolution of about 20 fim. As the detection 
of protons elastically scattered at angles down to few yurads requires a detector active area as 
close to the beam as ~ 1 mm, a novel "edgeless planar silicon" detector technology has been 
developed for TOTEM RPs in order to minimize an edge dead zone to only about 50 /im [I]. 

Each Tl telescope arm, covering the range 3.1 < \rj\ < 4.7, consists of five planes formed by 
six trapezoidal "Cathode Strip Chambers" (CSC) [lj (Fig. [6], center). These CSCs, with 10 mm 
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Figure 6: Left: Silicon detectors hosted in one pot. Center: One half-arm of the Tl telescope. 
Right: The installation of one half-arm of the T2 telescope. 

thick gas gap and a gas mixture of Ar/C0 2 /CF 4 (40%/50%/10%), provide three measurements 
of the charged particle coordinates with a spatial resolution of about 1 mm. The anode wires 
(pitch of 3 mm) will also provide level-1 trigger information; the cathode strips (pitch of 5 mm) 
are rotated by ± 60° with respect to the wires. 

The T2 telescope [5], based on "Gas Electron Multiplier" (GEM) technology [6], extends 
charged track reconstruction to the range 5.3 < |r/| < 6.5. Ten aligned detector planes, with a 
semicircular shape, are combined to form one of the half-arms located at ~ 13.5 m from IP5 
(Fig.0 right). This novel gas detector technology is ideal for the T2 telescope thanks to its good 
spatial resolution, excellent rate capability and good resistance to radiation. The T2 GEMs 
are characterized by a triple-GEM structure and a gas mixture of Ar/C0 2 (70%/30%) [5]. 
The read-out board has two separate layers with different patterns: 256x2 concentric circular 
strips (80 /im wide, pitch of 400 /mi), allow the track radial coordinate reconstruction with a 
resolution of about 100 /im; while a matrix of 24x65 pads (from 2x2 mm 2 to 7x7 mm 2 in size) 
provide level-1 trigger information and track azimuthal coordinate reconstruction. 

7 Physics Programme 

The physics goals of the TOTEM experiment are the measurement of the total pp cross section 
{ctot) with an ultimate precision of l-=-2%, the study of the nuclear elastic pp differential cross 
section (da e i/dt) over a wide range of \t\ (~ 10~ 3 < \t\ < 10 GeV 2 ) and the study of the inelastic 
interactions by measuring the cross section of soft diffractive processes and the forward charged 
particle flow. TOTEM can perform this programme operating in stand-alone mode, while the 
cooperation with CMS will make possible the study of hard diffraction, low-x physics, central 
exclusive diffractive production and the combination of particle multiplicity and energy flow in 
the forward region [2] . TOTEM measurements will allow to distinguish among different models 
of soft proton interactions, giving a deeper understanding of the proton structure. 

Fig. [7] (left), summarizing the existing measurements on pp and pp scattering, shows predic- 
tions for a tot from the COMPETE Collaboration based on fits according to different models [7j. 
The best fit predicts a to t = 111-5 ± 1.2l|j mb for the LHC energy (y/s = 14 TeV). The large 
uncertainties on available high energy data give a big error (90-=-130 mb), depending on the 
model used for the extrapolation. TOTEM aims to measure a to t with a precision down to 
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Figure 7: Left: Fits from the COMPETE Collaboration to available pp and pp scattering and 
cosmic ray data. Right: da e i/dt at the LHC as predicted by different models; t-acceptance 
ranges for different machine optics are also shown. 



l-j-2%, therefore allowing to discriminate among the different models. The measurement will 
be based on the "luminosity independent method" which, combining the optical theorem with 
the total rate, gives a to t and the machine luminosity (£) as a function of measurable rates: 

16tt dN el /dt\ t=0 = 1 + p 2 (N el + N inel ) 2 

° tot 1 + p 2 ' N el + N inel 16vr ' dN el /dt\ t=0 [ } 

where N e i and N ine i are respectively the elastic and inelastic rate and p is the ratio of real 
to imaginary part of the forward nuclear elastic scattering amplitude (given by theoretical 
predictions). Therefore, TOTEM will also provide an absolute measurement of C to be used 
for calibration purposes. The uncertainty on the extrapolation of dN^/dt to t = (optical 
point) depends on the acceptance for protons scattered at small \t\ values, hence at small 
angles. This requires a small beam angular divergence at the IP, which can be achieved in 
special runs with high (3* machine optics and typically low C. An approved optics with /3* = 
1540 m (and £ ~ 10 28 cm _2 s _1 ) will give a a tot (£) measurement at the level of l-i-2% (2%). 
Another approved optics with (3* = 90 m (and C ~ 10 30 cm _2 s _1 ), achievable without modifying 
the standard LHC injection optics, is expected to allow a preliminary o to t (£) measurement at 
the level of ~ 5% (7%) in the first period of the LHC running. The experimental systematic 
error for the measurement with j3* = 90 m will be dominated by the evaluation of dN e i/dt\ t= o, 
while with 0* = 1540 m it will be dominated by the uncertainty on the corrections to trigger 
losses in Single Diffraction events for masses below ~10GeV/c 2 [3J. Given the high rates 
involved, the statistical error on a tot will be negligible after few hours of data taking even at 
low C. The theoretical uncertainty related to the estimate of the p parameter is expected to 
give a contribution on the relative uncertainty of less than 1.2% (considering for instance the 
full error band on p extrapolation as derived in ref [7]). 

Fig. [7J (right) shows the distributions of da e i/dt at y/s = 14 TeV as predicted by different 
models in the whole |£|-range accessible by TOTEM according to the different LHC optics set- 
tings [3J. Several regions are found at increasing \t\: for \t\ < 6.5 x 10~ 4 GeV 2 (the Coulomb 
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region) da e i/dt ~ l/|t| 2 is dominated by photon exchange; for \t\ up to ~ 10~ 3 GeV 2 , the 
hadronic and Coulomb scattering interfere; for 10~ 3 < \t\ < 0.5 GeV 2 there is the hadronic 
region, e.g. described by "single-Pomeron exchange", characterized by an approximately ex- 
ponential fall (da e i/dt ~ e _B '*'); the diffractive structure of the proton is then expected for 
0.5 < \t\ < 1 GeV 2 ; for \t\ > 1 GeV 2 elastic collisions are described by pQCD, e.g. in terms of 
triple-gluon exchange with da e i/dt ~ \t\~ s - TOTEM will allow to discriminate among different 
models with a precise measurement of da e i/dt over all the accessible t-region, where da e i/dt 
spans over 11 orders of magnitude. In the hadronic region, important for the extrapolation of 
da e i/dt to t = 0, a fit on -B(|t|) is typically performed in the \t\ m i n < \t\ < 0.25 GeV 2 range, 
\t\min depending on the acceptance for protons scattered at small angles, hence on the beam 
optics. Fig. E] (left) shows the proton acceptance for the RPs at 220 m as a function of t for 
three different running scenarios, where the 50% acceptance is marked for each /3* value. For 
safety reasons, the minimum distance between RP detectors and the LHC beam is 10<7b eam ; by 
adding an extra 0.5 mm to take in account the distance from the edge of the sensitive detector 
area to the bottom of the RP window, this distance is about 1.3 mm for RP220 with /3* = 1540 
m, corresponding to a minimum angle of about 5 /zrad or equivalently 10~ 3 GeV 2 of squared 
4-momentum transfer. 

Diffractive (due to colour singlet exchange) and non-diffractive (due to colour exchange) 
inelastic interactions represent a big fraction (around 70 -j- 75 %) of a tot- Nevertheless many 
details of these processes, with close ties to proton structure and low-energy QCD, are still 
poorly understood. The majority of diffractive events exhibits intact ("leading") protons char- 
acterized by their t and fractional momentum loss £ = Ap/p. TOTEM will be able to measure 

t- and mass-distributions with acceptances depending on the beam optics. The diffractive 
proton acceptance for the RPs at 220 m has been computed for different optics as a function of 
£ and t (Fig. [HI right). For /3* = 1540 m or 90 m most of the protons are seen independently of 
their momentum loss. For low /3* values, elastic scattering is detectable only for \t\ > 2 GeV 2 , 
while diffracted protons are seen for £ > 2% for all \t\ values. 

When leading protons are detected on both sides of IP5, as in the case of Double Pomeron 
Exchange (DPE), the central mass M acceptance is shown in Fig. [HI left. By combining data 
from low, intermediate and high (3* runs, the differential cross-section as a function of M can 
be measured with good precision over the full mass range. 

8 Early Physics Programme 

LHC will soon restart with a beginning center of mass energy of 7 TeV and a low 0* optics. 
We are studying the TOTEM performance under these conditions, however we are also asking 
for /3* =90 m for short runs in 2010. The acceptance for protons at the 220 m RP stations is 
shown in Fig. [9] (right); they are mainly seen for a momentum loss £ > 2%, with a resolution 
on £ slightly less than 10 -2 . These results are for a low (3* optics, they depend very strongly on 
the beam optics and little on the energy, so that they are similar to those obtained for 14 TeV. 
For DPE, the predicted central diffractive mass distribution would start at around 250 GeV, a 
significant difference from the very low mass of few GeV reachable with a high /3* optics (see 
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Figure 8: Left: Proton acceptance for the RPs at 220 m as a function of t for three different 
running scenarios, where the 50% acceptance is marked for each /3* value. Right: Diffractive 
proton acceptance for the RPs at 220 m for different optics as a function of £ and t. 




Figure 9: Left: The DPE central mass M distribution for three (3* values. Right: Proton 
acceptance for early runs at the 220 m RP stations. 

Fig. [91 left). Therefore the very first measurements with a low 0* will consist in: 1) by using 
the horizontal RPs, in Single Diffractive and DPE for high masses of the diffractive system; 2) 
by using the vertical RPs for the measurement of the elastic scattering over a range of 1< \t\ < 
10 GeV 2 , which would not allow a measurement of a to t but, however, it will already allow to 
discriminate among different models for the prediction of da e i/dt. 

The Tl and T2 telescopes will provide a measurement of charged particle multiplicity for 
different processes as well as the identification and measurement of rapidity gaps. Primary 
cosmic rays in the PeV (10 15 eV) energy range and above are a challenging issue in astrophysics. 
The LHC center of mass energy corresponds to a 100 PeV energy for a fixed target collision in the 
air, at the same time providing a high event rate relative to the very low rate of cosmic particles 
in this energy domain. A primary cosmic ray entering the upper atmosphere experiences 
a nuclear interaction, with the production of nuclear fragments and ir mesons, starting an 
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air shower with hadronic, electromagnetic and muon components. The real challenge is to 
determine the nature of the primary interaction and the energy and composition of the incident 
particle from the measurement of the shower. Several high energy hadronic interaction models 
are nowadays available, which predict energy flow, multiplicity and other quantities of such 
showers. There are large differences between the predictions of currently available models, with 
significant inconsistencies in the forward region. Among the several quantities that can be 
measured by TOTEM and compared with model predictions, the charged particle multiplicity 
in the Tl and T2 acceptance region shows significant differences in the predictions obtained 
by different available hadronic interaction models for cosmic ray showers, once the events are 
passed through the simulation of Tl and T2. Therefore this measurement can already be used 
in early runs to validate/tune the generators [2]. 

9 Summary and Conclusions 

The TOTEM experiment will be ready for data taking at the LHC restart. Running under all 
beam conditions, it will perform an important and exciting physics programme involving the 
measurement of a to t and da e i/dt in pp interactions as well as studies on diffractive processes 
and on forward charged particle production. Special high (3* runs will be required in order to 
measure a to t at the level of ~ 5 % (early measurement with 0* = 90 m) and ~ 1 4- 2 % (/?* = 
1540 m). da e i/dt will be studied in the range ~ 10~ 3 < \t\ < 10 GeV 2 allowing to distinguish 
among several theoretical predictions. A common physics programme with CMS on soft and 
hard diffraction as well as on forward particle flow studies is also foreseen in a later stage. 
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Abstract 

The ATLAS forward detector system is presented. Following this, the forward physics 
measurements that are expected to be carried out with early ATLAS data are introduced 
and the relevant trigger and analysis strategies discussed. 

1 The ATLAS forward detector system 

ATLAS is a general multipurpose detector P designed to be able to measure a wide variety of 
physics processes at the LHC. ATLAS has a standard layout of sub-detectors; an Inner Detector 
(\r)\ < 2.5) for tracking purposes, electromagnetic calorimeters (\rj\ < 3.2) for measuring the 
energy of electrons and photons, hadronic calorimeters (\r]\ < 4.9) to measure the energy of 
mesons and baryons, and a muon spectrometer (\rj\ < 2.7). In addition, there are a number of 
other sub-detectors designed to measure (forward) particle production; these are the MBTS, 
LUCID, ZDC and ALFA. 

The Minimum Bias Trigger Scintillators (MBTS) are located at \z\ = 3.6 m and cover 
2.1 < 1 7] | < 3.8 [2]. The MBTS is designed to trigger so-called minimum bias events, identifying 
particle production from inelastic collisions. Each side of the MBTS (positive and negative) are 
divided into 16 segments: two 77-rings and eight 0-sections. The multiplicity of hit segments 
on each side will be available for the Level 1 (LI) trigger. The MBTS detector will only be 
available during the low luminosity phase of LHC operation, as it was not designed to have 
high radiation tolerance. 

The LUCID detectors [3] are located 17 m from the interaction point, one on each side of 
ATLAS, and provide coverage of 5.6 < < 6.0 for charged particles. Each LUCID detector 
is a symmetric array of polished aluminium tubes that surround the beam-pipe. Each tube is 
15 mm in diameter and filled with C4 F10 gas, which results in Cerenkov emission from charged 
particles crossing the tube. The Cerenkov light is read out by photo-multiplier tubes. The 
initial LUCID detector, available for early data taking, has LI trigger capabilities but only 
limited azimuthal coverage. An upgrade to provide full azimuthal coverage is under study. 

The Zero Degree Calorimeter (ZDC) [4] is located 140 m from the interaction point in the 
TAN region (target absorber for neutrals), where the single beam-pipe splits into two, and 
provides coverage of \r)\ > 8.3 for neutral particles. The ZDC consists of one electromagnetic 
and three hadronic tungsten/quartz calorimeters. Vertical quartz strips provide the energy 
measurements and horizontal quartz rods are used for coordinate readout. At LHC startup, 
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when there are few bunches in the beam, the electromagnetic calorimeter is not installed and 
the space it would occupy is used by the LHCf experiment. After initial running, LHCf will be 
removed and the full ZDC installed. 

The ALFA roman pot (RP) spectrometers are located 240 m from the interaction point [5]. 
Unlike other detectors, the RP spectrometers are not fixed relative to the beam. At injection, 
the ALFA detectors are in a withdrawn position far from the beam. After the protons have 
been injected and the beam has been stabilized, the detectors are moved to within 1.5 mm of 
the beam. Elastic and diffractive protons that have been deflected outside the beam envelope 
pass through arrays of scintillating fibre trackers (20x64 fibres in each array), which measure 
the position of the protons with respect to the beam. This allows the momentum and angle 
of each proton to be reconstructed by using the known LHC lattice. ALFA will only be used 
during special LHC runs at low luminosities with high /3* optics. 



2 Soft diffraction 

Single diffraction is a low t-process in which a colour singlet (i.e. pomeron) is exchanged between 
the two protons and one of the protons breaks up into a dissociative system. There is a large 
rapidity gap between the outgoing proton and the dissociative system as a colour singlet object 
was exchanged. The cross section for soft single diffraction (SD) is expected to be of the order 
10 mb at the LHC. There are two approaches that will be used to measure soft-SD at ATLAS. 

The first approach will identify the dissociated system using the inner detector, calorimeters, 
LUCID and possibly the ZDC. The variable of interest is the fractional longitudinal momentum 
loss, £, given by 

where Mx is the invariant mass of the dissociative system and s is the centre-of-mass energy of 
the pp collision. It follows that the dissociated system for events with low-£ will be contained 
only in the forward detectors, whereas high-£ events will have activity in many areas of the 
central detector as well. The MBTS, LUCID and the ZDC will be required to trigger events 
across the full kinematic range. It has been estimated that a sample of one million events can 
be collected in two weeks given a luminosity of 10 31 cm" 2 s . 

The second approach is to tag the outgoing proton and measure £ directly using 

* = 1 - FT ( 2 ) 

where p z and p' z are the longitudinal momenta of the incoming and outgoing protons respectively. 
This approach requires the ALFA RP detectors, which can only be used in special runs with 
high-/3* optics at a luminosity of 10 27 cm~ 2 s -1 . ALFA will be able to measure the fractional 
momentum loss, with an accuracy between 8% (for £ ~ 0.01) and 2% (for £ ~ 0.1). The trigger 
will be provided by ALFA and it is expected that 1.2-1.8 million events will be retained with 
just 100 hrs of data taking [5]. The analysis will also make use of LUCID and the ZDC, which 
are used to tag the dissociative system to separate the events from elastic scattering. 
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In addition to soft single diffraction, it is expected that ATLAS will also be able to measure 
soft double diffraction, in which both protons form dissociative systems. 

3 Diffractive di-jet production 

ATLAS will be able to measure single diffractive di-jet production and di-jet production via 
double pomeron exchange. These processes, shown in Fig. [1] (a) and (b), will allow the study 
of diffractive parton density functions (dPDFs) and factorization breaking in diffractive events. 
Factorization breaking is the observation that the dPDFs obtained at HERA do not predict the 
correct cross section for diffractive events at hadron colliders [6] and is attributed to secondary 
scattering between spectator partons in the protons causing the proton to break up and the 
rapidity gap to be destroyed. Factorization breaking will be studied by measuring the ratio of 
single diffractive to non-diffractive di-jet events, R(SD/ND), and the ratio of double-pomeron 
exchange to single diffractive di-jet events, R(DPE/SD). 




(a) (b) (c) 

Figure 1: Single diffractive di-jet production (a), di-jet production via double pomeron exchange 
(b) and central exclusive di-jet production (c). In (a) and (b), the 'zig-zag' line represents colour 
singlet (pomeron) exchange and there will be rapidity gap between the intact proton and the 
pomeron remnants. In (c), the di-jets are produced exclusively, implying no hadronic activity 
outside of the di-jet system. 

It is expected that a few thousand SD di-jet events, with jet transverse energy, > 20 GeV, 
will be available for analysis given 100 pb _1 of data. The number of events is restricted by the 
large prescale applied in the Level 1 (LI) trigger for low transverse energy jets and the current 
focus is to develop new trigger strategies capable of retaining more events. The new triggers 
may include information from LUCID and/or the MBTS (discussed futher in Sec. H|). 

4 Central exclusive di-jet production 

Central exclusive di-jet production (CEP) is the process pp — >■ p + jj+p, where the '+' denotes 
a large rapidity gap from the outgoing protons, and is shown in Fig. [1] (c). CEP has received 
a great deal of attention in recent years due to the possibility of tagging the outgoing protons, 
using new forward proton detectors, in order to measure properties of the Higgs boson [TJ. 
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Although the CDF measurements of exclusive di-jets, photons and charmonium are in good 
agreement with the theoretical predications [8], there remains a factor of three uncertainty 
in the calculation of the Higgs cross section at 14 TeV. Exclusive di-jet production offers the 
opportunity to constrain this uncertainty at the LHC with early ATLAS data. 

The analysis strategy is to define the exclusivity of an event using the di-jet mass fraction, 

«* = P~ (3) 

where My is the invariant mass of the di-jets and M ca \ in the mass of all energy deposits in 
the calorimeter. Typically, an exclusive event will have ~ 1 and inclusive /diffract ive events 
will have R-^ <C 1. However, the tail of the background distribution extends up to large 
values and the extraction of the exclusive signal will require a very good understanding of the 
background distributions. 

Due to the large QCD cross sections, the standard low-.&r jet triggers will be heavily 
prescaled and will provide insufficient statistics to perform CEP measurements. Consequently, 
a new trigger was developed that, in addition to a jet with E T > 18 GeV, required an absence 
of hit segments on at least one side of the MBTS. Simulations show that this trigger has a 60% 
efficiency for signal events (with respect to the jet trigger alone) and provides a 10 4 rejection 
factor for standard QCD events. It is expected that a few hundred exclusive di-jet events will 
survive the final selection criteria for every 10 pb _1 of data. 

In addition to exclusive studies, this trigger strategy can also be used to retain single 
diffractive di-jet events. Although the overall trigger efficiency for single diffractive events is 
approximately 10 -2 (for this trigger), the majority of the events that are removed are at high-£. 
Thus, the trigger may allow a sample of low-£ SD di-jet events to be analyzed, which could 
then be used to constrain the dPDFs in an interesting kinematic region. 



5 Forward jets 

There are a number of signatures involving widely separated jets that can be measured with 
early ATLAS data. The processes of interest are 2 — > 2 partonic scatters mediated by t-channel 
colour octet or colour singlet exchange. 

An inclusive measurement of interest is the azimuthal de-correlation of di-jets, the size of 
which is dependent on the pseudo-rapidity separation of the jets . The theoretical predictions 
of the azimuthal de-correlation is dependent on the framework used to make the prediction, 
namely fixed order (NLO), parton shower and BFKL. It is expected that LHC measurements 
will be sensitive to BFKL effects [9]. 

It is also useful to split the di-jet samples into gap and non-gap components. This is 
achieved by vetoing on radiation between the jets. Experimentally, we introduce a veto-scale 
which defines the maximum transverse energy that can be deposited between the jets in a gap 
event. It is then useful to study the fraction of gap-events as a function of the pseudo-rapidity 
separation of the jets, Arj. For large separation, At] ~ 6, the gap events are predominantly due 
to colour singlet exchange. A prediction of BFKL [10] is that the fraction of events with little 
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activity between the jets should rise with the separation of the jets, Arj, and was extensively 
searched for at the Tevatron and HERA [TT]. The rise of the gap- fraction was not observed at 
the Tevatron, for example, because the centre-of-mass energy was too small; it was shown in 
[12] that the rapidly falling PDFs at high x tempered the rise and meant that a large enough 
sample of events with large jet separations could not be obtained (because A77 = ln(s/t)). 
An improved measurement should be possible at the LHC due to the increased centre-of-mass 
energy. In principle, ATLAS should be able to measure the gap-fraction up to Arj ~ 9, 9.5 with 
approximately 10 pb _1 of data. The events are retained for analysis by using a forward di-jet 
(|r/ jet | > 3.2) trigger with E T > 18 GeV. 

Finally, for smaller pseudo-rapidity separations (Arj ~ 4), it is interesting to study the 
fraction of gap events as a function of ln(Q/ Qo), where Q is the transverse energy of the jets. In 
this region, the gap events also arise from single gluon (colour octet) exchange. This observable 
is sensitive to QCD effects such as wide angle soft gluon radiation [131 EH]- Preliminary studies 
indicate that these measurements should be achievable with less than 100 pb _1 of data using 
standard jet triggers. 
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Abstract 

We describe a new framework for calculating multiple hard jet emission in Higgs 
boson production via gluon-gluon fusion. It is based on the well-known FKL factorisation 
formula. This formally applies only in a certain high energy limit, but we describe how to 
modify the description so as to work well at LHC energies. Approximate matrix elements 
thus obtained are validated by comparing with known tree level results, and the framework 
is implemented in a Monte Carlo event generator. 

1 Introduction 

Scattering events at the LHC will be dominated by multijet final states, ultimately arising 
from QCD radiation. The scattering amplitudes for such events involve Feynman diagrams 
with many external quark and gluon legs, which cannot be feasibly calculated using present 
computing power. Instead, approximate methods are usually used to estimate multijet final 
states. A common practice is to interface tree level matrix elements with parton shower al- 
gorithms. These simulate the effect of QCD radiation on a given hard matrix element, where 
the extra partons are approximated well in the soft and collinear limit i.e. when the emitted 
radiation has a low transverse momentum with respect to the external legs of the hard in- 
teraction. Such an approach, whilst undoubtedly useful, is known to be inadequate in many 
processes. The question then naturally arises of whether it is possible to estimate scattering 
amplitudes with many outgoing hard partons directly (i.e. with no restriction on the transverse 
momentum). 

One such process in which multiple hard jet production is relevant is Higgs boson production. 
There are two main production modes at the LHC, known as W boson fusion (WBF) and gluon- 
gluon fusion (GGF) and depicted in figure [H Both of these can be used as a discovery channel, 




WBF GGF 



Figure 1: The two main Higgs production modes at the LHC. 

and each forms a background to the other. In WBF there is no colour exchange in the t-channel, 
and any jet activity is mostly confined to the forward and backward detector regions. In GGF 
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however, in which the Higgs couples to gluons via a top quark loop, there is a colour octet 
exchange. One thus expects jet activity over the complete range in rapidity, and the nature 
of the radiation (e.g. whether it is hard or soft) has a significant impact on what one can 
measure about the Higgs boson, as we will see. Thus, a means for accurately estimating hard 
jet radiation in GGF is extremely well-motivated. 



2 FKL factorisation 

Our aim is to estimate scattering amplitudes with many external hard parton legs, and our start- 



ing point is the FKL factorisation formula of Fadin et al.(1975)Fadin, Kuraev, and Lipatov 



This states that in the limit of multi-Regge kinematics (MRK) in which the rapidity gap between 
final state particles is large (i.e. large centre of mass energy but fixed momentum transfer), the 
amplitude for pure multiparton scattering (from initial state partons a and 0) is dominated by 
the process 

a + f3 ->■ a + f3 + ng, (1) 

for a given number (n + 2) of final state partons. Furthermore, the sum of these diagrams gives 
a factorised expression for the scattering amplitude, which (including also a Higgs boson in the 
final state) has the form 



(m+i . 
i=l qi 

' ni \ 

(ft-i>?i) J C H (q ni+ i,q ni . 

n+l 

II — e ^v[a{q 2 i ){y l -i-y l )} 

=ni+2 qi 

J] C^q^qA. (2) 

=ni+2 / 

In this formula, are the 4-momenta of the virtual gluons, Ch is an effective vertex coupling 
the Higgs to gluons |Del Duca et al.(2003)Del Duca, Kilgore, Oleari, Schmidt, and Zeppenfeld 



C M the Lipatov effective vertex for gluon emission, and a(q 2 ) a function which includes leading 



virtual corrections in the MRK limit (see Andersen and White(2008) , Andersen et al. (2009) Andersen, Del Du 



for more details). The form of this formula is shown schematically in figure [2j 

The chief advantages of this formula are firstly that there is no restriction on the ordering of 
final state transverse momenta. Secondly, the factorised form means that scattering amplitudes 
with any number of partons can be efficiently calculated, thus this is indeed a good starting 
point for trying to describe hard final state radiation. Furthermore, the virtual corrections 
built into the Reggeised gluon propagators cancel the collinear singularities arising from real 
gluon emission. However, a significant disadvantage of eq. (T5]) is that it formally applies only 
in the MRK limit, which is not well-approximated by the phase space probed at the LHC or 
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Figure 2: Diagrams contributing to the GGF scattering amplitude in the MRK limit, where the 
gluons couple with Lipatov vertices, and the virtual gluons are Reggeised (i.e. have resummed 
leading virtual corrections). 

Tevatron. The solution is to modify the application of the FKL formula so as to build in known 
features of perturbation theory. 



3 Modified approach 

Our approach to calculating Higgs boson plus multiparton scattering amplitudes can be sum- 
marised as follows. We define the scattering amplitude with (n + 2) final state partons by 
eq. ([2]) implemented according to the following prescription: 

1. Impose 4- momentum conservation at the emission vertices. 

2. Use full virtual 4-momenta instead of transverse components only. 

3. Require gauge invariance of the Lipatov vertex (k ■ C=0) over all of phase space. 

The first requirement is an obvious piece of physics, but is formally subleading in the cen- 
tre of mass energy so is not included in most previous analytic studies of eq. d2J). It im- 
plements corrections to the matrix element, as well as the phase space of the emitted glu- 
ons. The replacement of virtual 4-momenta by their transverse components (as is usually 



done in the BFKL implementation of eq. fl2]) Fadin et al.(1975)Fadin, Kuraev, and Lipatov 



Kuraev et al.(1976)Kuraev, Lipatov, and Fadin Kuraev et al.(1977)Kuraev, Lipatov, and Fadin] ) 



is allowable in the MRK limit, but leads to significant differences away from this corner of phase 
space. Using the full 4-momenta amounts to keeping known singularities of the scattering ampli- 
tude in the same place outside the MRK limit, rather than shifting them. The final requirement 
has the consequence that the squared Lipatov vertex is positive, thus enforcing local positivity of 



the matrix element (see Andersen and White(2008), Andersen et al. (2009) Andersen, Del Duca, and White 



for details). The gauge invariance constraint overlaps with, but is not the same as, the kinematic 



constraint of Catani et al.(1990)Catani, Fiorani, and Marchesini . 
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Figure 3: Comparison of the cross-sections for Higgs boson production via GGF in association 
with two and three partons obtained using the modified FKL approach, and the full tree level 
results from MadGraph. Uncertainty corresponds to scale variation by a factor of two. 



Each of the above requirements goes beyond any logarithmic order in the high energy 
expansion, and the resulting matrix elements are in a sense closer to fixed order perturbation 
theory than to the original FKL description, as we will see. 

In order to judge whether the approximate matrix elements work well, one may compare 
them with the known tree level results for Higgs boson production in association with two and 
three partons, which we have obtained from MadGraph Alwall et al.(2007)] . This comparison is 
made in figure [3]for WBF-like signal cuts (see Andersen et al. (2009) Andersen, Del Duca, and White] ). 
These have been applied inclusively i.e. at least one pair of jets is required to pass the cuts. 
One sees that the modified FKL results are well within the scale variation of the tree level 
results. Alternatively, the results obtained using a traditional BFKL implementation (but also 
including energy and momentum conservation) are shown in figure HI One sees that this lies 
outside the scale-variation associated with the tree level results. Without 4-momentum conser- 



vation, the situation is far worse Andersen et al. (2009) Andersen, Del Duca, and White , and 
indeed worsens at higher orders due to the ability to create multiple gluons at no cost in energy. 



Kinematic distributions as well as total cross-sections are well approximated using the modi- 
fied FKL approach. As examples, we show the transverse momentum and rapidity distributions 
of the Higgs boson in figures |5] and El for the case of Higgs boson production in association 
with three partons. One sees that the results obtained from the modified FKL approach agree 
closely with the full tree level results, and are certainly well within the scale variation uncer- 
tainty. Another example is given in figure which shows the distribution of azimuthal angle 
between the tagging jets. The shape of this distribution is intimately related to the CP nature 



of the Higgs boson Hankele et al.(2006)Hankele, Klamke, and Zeppenfeld . Again, very good 
agreement is observed between the modified FKL description and the full tree level result. 

Having validated the approximation by comparing to fixed order results, one may then 
use the modified FKL formalism to estimate multiparton matrix elements with any number 
of final state partons. We have implemented the framework in a Monte Carlo event gener- 
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Figure 4: Comparison of the cross-sections for Higgs boson production via GGF in association 
with two and three partons obtained using a traditional BFKL approach (with 4-momentum 
conservation implemented), and the full tree level results from MadGraph. Uncertainty corre- 
sponds to scale variation by a factor of two. 




Figure 5: Distribution of the Higgs boson transverse momentum for the hjjj final state. 




Figure 6: Distribution of the Higgs boson rapidity, for the hjjj final state. 
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Figure 7: Distribution of the azimuthal angle between the tagging jets, for the hjjj final state. 

ator for Higgs boson plus multijet production. The tree level matrix elements for hjj and 
hjjj production are also included, with a matching procedure to avoid any double count- 
ing 



Andersen et al. (2009) Andersen, Del Duca, and White , and code is available at 



http : / / andersen . web . cern . ch/ andersen/M JEV/ 
index.html 

The known tree level results (i.e. for two and three final state partons) have also been included 



using a suitable matching procedure to avoid double counting Andersen and White (2008) 



Andersen et al. (2009) Andersen, Del Duca, and White , and here we show example results. In 



figure Ufa) we show the resummed distribution in the number of hard jets, obtained using 
inclusively applied WBF-like signal cuts, as above. In figure Mb) we show the same distribution, 
but where the hardest jets are explicitly required to pass the cuts. One sees in both cases that 
there are a significant number of events with more than three hard jets. Furthermore, when the 
hardest jets are used as the tagging jets, the impact of higher order jet radiation is reduced. 

The azimuthal angle <j)j a j b between the tagging jets is depicted in figure |9] (i.e. the fully re- 
summed equivalent of figure [7]). The choice of cuts has a significant impact on how much decorre- 
lation is observed, and the hardest jet cuts lead to a stronger observed correlation than the inclu- 
sive cuts (consistent with figured]). Further phenomenological investigation is necessary in order 
to determine whether the CP nature of the Higgs boson can be fully determined i.e. whether 
a strong enough correlation can be preserved without cutting into the cross-section too much. A 



useful quantity to consider is the parameter jHankele et al.(2006)Hankele, Klamke, and Zeppenfeld 



A = V{<t>0«h < ^Z 4 ) - < hah < 3 V 4 ) + <<t>iai b > 37T/4) 

* ^ jajb < tt/4) + a(vr/4 < cj> jajb < 3vr/4) + a{<f> jajb > 3tt/4) ' 1 ) 

where a is the fully integrated cross-section. This parameter is positive in the case of a 
scalar Higgs boson, and negative for a pseudo-scalar. Thus, a measurement of the correla- 
tion parameter A$ allows one to infer whether or not a discovered Higgs boson is of Stan- 
dard Model character or not. Values of obtained in different calculations are shown 
in table [31 The results labelled LO correspond to the tree level calculations for hjj and 
hjjj production, obtained using MadGraph. Results are shown for WBF-like signal cuts 
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Figure 8: Resummed distribution in the number of hard jets i.e. those jets satisfying p t > 
40GeV, in the cases where (a) the WBF-like cuts are applied inclusively; (b) the hardest jets 
are required to pass the cuts. Jets have been clustered with a k t algorithm, and uncertainty 
corresponds to scale variation by a factor of two. 



Inclusive cuts 




Hardest cuts 




LO 2-jet 


0.456 


LO 2-jet 


0.456 


Resummed, = 2-jet 


0.444 


Resummed, = 2-jet 


0.436 


LO 3-jet 


0.203 


LO 3-jet 


0.374 


Resummed 


0.123 


Resummed 


0.372 



Table 1: Values of the correlation parameter of eq. ([3]) using tree-level matrix elements (first 
and third rows), and the resummed approach described here. Results are shown for the two 
choices of cuts discussed in the text. 



(see jAndersen et al. (2009) Andersen, Del Duca, and White] ) , which are applied both inclusively 
(left-hand columns) and in the case where the hardest jets are used as the tagging jets (right- 
hand columns). One sees that the degree of correlation observed is highly dependent on the 
cuts used, and that much more decorrelation is observed for the inclusive cuts. Furthermore, 
the emission of one extra hard parton already leads to significant decorrelation, which is then 
exacerbated by the resummation (albeit less severely when the hardest jets are used as the 
tagging jets). Finally, one recovers a high degree of correlation by requiring that only two hard 
jets are present, even when the full resummation is included. 



4 Conclusions 

We have described a framework for estimating scattering amplitudes with many hard final 
state partons, and applied this to Higgs production via gluon-gluon fusion at the LHC. The 
approach is based on the FKL factorisation formalism, which is then modified to incorporate 
known features of perturbation theory. The results have been validated by comparing with 
known tree level results for the two and three parton final states, and a Monte Carlo generator 
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Figure 9: Resummed distribution of the azimuthal angle between the tagging jets, in the case 
where (a) the WBF-like signal cuts are applied inclusively; (b) the hardest jets are required to 
pass the cuts. 

for GGF is available. 

The approach described here is not limited to Higgs boson production, but is also applicable 
to other processes such as W production in association with jets, and pure QCD jet production. 
Furthermore, the underlying FKL factorisation formula can be improved, although the leading 
order description is already working rather well. 

It would also be interesting to interface the modified FKL approach with a parton shower 
algorithm. The former technique is complementary to the latter in a well-defined sense, and 
should estimate the final state jet topology well rather than the jet substructure (which is 
approximated well by the parton shower). 
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Studying the Underlying Event at the Tevatron 



Deepak Kar, On behalf of the CDF Collaboration 
Abstract 

CDF Run II data for the underlying event associated with Drell-Yan lepton pair pro- 
duction are studied as a function of the lepton-pair transverse momentum. The data are 
compared with a previous analysis on the behavior of the underlying event in high trans- 
verse momentum jet production and also with several other QCD Monte-Carlo models. 
The goal is to provide data that can be used to tune and improve the QCD Monte-Carlo 
models of the underlying event, which is especially important now for the startup of the 
Large Hadron Collider. 

5 Introduction: the Underlying Event 

In order to find 'new' physics at a hadron-hadron collider it is essential to have Monte-Carlo 
models that simulate accurately the 'ordinary' QCD hard-scattering events. To do this one 
must not only have a good model of the hard scattering part of the process, but also of the 
underlying event. 

A typical 2-to-2 hard scattering event is a proton-antiproton collision at the hadron colliders 
as shown in the Fig. 1, all happening inside the radius of a proton. In addition to the two hard 
scattered outgoing partons, which fragment into jets - there is initial and final state radiation 
(caused by bremsstrahlung and gluon emission), multiple parton interaction (additional 2-to-2 
scattering within the same event), 'beam beam remnants' (particles that come from the breakup 
of the proton and antiproton, from the partons not participating in the primary hard scatter). 
We define the 'underlying event' pQ as everything except the hard scattered components, which 
includes the 'beam-beam remnants' (or the BBR) plus the multiple parton interaction (or the 
MPI). However, it is not possible on an event-by-event basis to be certain which particles came 
from the underlying event and, which particles originated from the hard scattering. The 'un- 
derlying event' (i.e. BBR plus MPI) is an unavoidable background to most collider observables. 
For example, at the Tevatron both the inclusive jet cross section and the b-jet cross section, as 
well as isolation cuts and the measurement of missing energy depend sensitively on the under- 
lying event. A good understanding of it will lead to more precise measurements at the Tevatron 
and the LHC. 

Experimentally it is possible to take advantage of the topological structure of hadron-hadron 
collisions to study the underlying event. The direction of the Z boson is used to isolate regions 
of f] — <p space that are sensitive to the underlying event. The angle A0 = — (f>z is the relative 
azimuthal angle between charged particles coming from the underlying event and the direction 
of Z boson, as in Fig. 2 (left). We split the central region defined between \r]\ < 1 as follows, 

• \Acj)\ < 60° as the toward region. 
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Figure 10: A typical 2-2 hard scattering process 




Figure 11: On the left, dividing the central region, with relative to the Z boson direction and 
on the right, TransMAX and transMIN regions. 

• 60° < \A(f)\ < 120° as the transverse region. And, 

• \A(j)\ > 120° as the away region. 

As illustrated in Fig. 2 (right), we define MAX and MIN transverse regions which help to 
separate the hard component (initial and final-state radiation) from the beam-beam remnant 
component. MAX (MIN) refer to the transverse region containing largest (smallest) number 
of charged particles or to the region containing the largest (smallest) scalar p-psum of charged 
particles, on an event by event basis. One expects that the transMAX region will pick up the 
hardest initial or final-state radiation while both the transMAX and transMIN regions should 
receive beam-beam remnant contributions. Hence one expects the transMIN region to be more 
sensitive to the beam-beam remnant component of the underlying event, while the transMAX 
minus the transMIN (i.e., transDIF) is very sensitive to hard initial and final-state radiation. 
This idea, was first suggested by Bryan Webber and Pino Marchesini [2], and implemented in 
a paper by Jon Pumplin [3]. 

For hard scattered jets the transverse regions are most sensitive to underlying events, since 
they are perpendicular to the plane of 2-to-2 hard scattering. For Drell-Yan its easy to identify 
and remove leptons (since they are the colorless components) from the transverse and toward 
(which can not be done for dijet events, as the leading jet is itself in toward region) regions and 
use them to study the underlying event. So we can see that Drell-Yan events are a clean probe 
of the underlying events. 
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6 Comparing data with QCD Monte Carlo Models 



6.1 The underlying event as a function of lepton pair px 

We looked at the charged particles in the range pt > 0.5 GeV/c and \rj\ < 1, at the region of 
Z-boson, defined as 70 GeV/c 2 < Mu < 110 GeV/c 2 , in the 'toward', 'away' and 'transverse' 
regions, as defined in Fig. 2. The electron and muon selections are based on the standard CDF 
high pt electron and muon selection criteria. The lepton pairs are formed by oppositely charged 
leptons, with the requirement that both the leptons came from the same primary collision. We 
use a time of flight [I] cosmic filter to eliminate cosmic muons, as the time difference between 
the muons recorded in the upper and lower half of the detector is expected to be very small 
for muons not coming from cosmic rays. The Z — > e + e~ and Z — » data sample contains 

backgrounds mainly from QCD jets and W+jets. Studies [5] have shown that these backgrounds 
are negligible at the region of Z boson. We use the ratio of the generator level Monte Carlo result 
[B] and the detector level Monte Carlo result as our correction factor for correcting the data 
back to the particle level. We analyzed data corresponding to the luminosity of approximately 
2.7 fb-\ 

The underlying event observables are found to be reasonably flat with the increasing lepton 
pair transverse momentum in the transverse and toward regions, but goes up in the away region 
to balance the lepton pairs. 

In Fig. 3, we looked at the two observables corresponding to the underlying event, the 
number of charged particle density and the charged transverse momentum sum density in the 
transverse region, compared with pythia [7] tunes A and AW [8], herwig [H] without MPI 
and a previous CDF analysis on leading jet underlying event results. We overlay the results 
for all the regions in Fig. 4. In Fig. 5 and 6, we look at the same observables in transMAX, 
transMIN, and transDIF regions. Overall, we can see that pythia tune AW does a good job 
of reproducing the data, herwig (without MPI) does not produce enough activity in the 
transverse region for either process. There is no final-state radiation in Z-boson production so 
that the lack of MPI becomes more evident, herwig (with jimmy [10J MPI) agrees with tune 
AW for the scalar pT sum density in the toward and transMIN regions. However, it produces 
too much charged particle density in these regions, herwig (with jimmy MPI) fits the pt 
sum density, but it does so by producing too many charged particles {i.e. it has too soft of a 
Pt spectrum in these regions). This can be seen in Fig. 7 which shows the data for Z-boson 
events on the average charged particle pt and the average maximum charged particle pt, for 
the transverse region compared with the QCD Monte Carlo models. So the pt distributions in 
the transverse region are too soft, resulting in an average pt and average pt maximum that are 
too small. Comparing HERWIG (without MPI) with HERWIG (with JIMMY MPI) clearly shows 
the importance of MPI in these regions. 

We also compared them with leading jet underlying event results [11 J and observed rea- 
sonably close agreement - which may indicate the universality of underlying event modeling. 
However, at low pt, we see a difference. If the leading jet has no transverse momentum then 
there are no charged particles, we just get min-bias events. There are a lot of low transverse 
momentum jets and for pt (jet#l) < 30 GeV/c the leading jet is not always the jet resulting 
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Figure 12: Drell-Yan underlying event plots, charged particle multiplicity at the top and the 
charged px sum at the bottom. 
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Figure 13: Overlaying all three regions, charged particle multiplicity at the top and the charged 
Pt sum at the bottom. 

from the hard 2-to-2 scattering. This produces is a bump in the transverse density at low px ■ 



6.2 Correlation studies 

The rate of change of < pr > versus charged multiplicity is a measure of the amount of hard 
versus soft processes contributing to collisions and it is sensitive to the modeling of the multiple 
parton interactions [12]. This variable is the most poorly reproduced variable by the available 
Monte-Carlo generators. If only the soft beam-beam remnants contributed to min-bias collisions 
then < pt > would not depend on charged multiplicity. If one has two processes contributing, 
one soft (beam-beam remnants) and one hard (hard 2-to-2 parton-parton scattering), then 
demanding large multiplicity would preferentially select the hard process and lead to a high 
< Pt >■ However, we see that with only these two processes < pt > increases much too 
rapidly as a function of multiplicity. Multiple-parton interactions provides another mechanism 
for producing large multiplicities that are harder than the beam-beam remnants, but not as 
hard as the primary 2-to-2 hard scattering. 

Fig. 8 (top) shows the data corrected to the particle level on the average pt of charged 
particles versus the multiplicity for charged particles with p T > 0.5 GeV/c and \r)\ < 1 for 
Z-boson events from this analysis. HERWIG (without MPI) predicts the < > to rise too 
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Figure 14: TransMAX and transMIN regions, charged particle multiplicity at the top and the 
charged pt sum at the bottom. 
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Figure 15: TransDIF regions, charged particle multiplicity at the top and the charged px sum 
at the bottom. 
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Figure 16: Drell-Yan underlying event plots, charged particle < pt > at the top and the charged 
Pt maximum at the bottom. 
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Figure 17: Charged multiplicity against charged transverse momentum average correlation 
plots. 

rapidly as the multiplicity increases. For HERWIG (without MPI) large multiplicities come from 
events with a high pt Z-boson and hence a large pt 'away-side' jet. This can be seen clearly 
in Fig. 8 (middle) which shows the average pt of the Z-boson versus the charged multiplicity. 
Without MPI the only way of getting large multiplicity is with high pr(Z) events. For the 
models with MPI one can get large multiplicity either from high pr(Z) events or from MPI and 
hence < Pt{Z) > does not rise as sharply with multiplicity in accord with the data, pythia 
tune AW describes the Z-boson data fairly well. Fig. 8 (bottom) shows the data corrected 
to the particle level on the average pt of charged particles versus the multiplicity for charged 
particles with p T > 0.5 GeV/c and \rj\ < 1 for Z-boson events in which pr{Z) < 10 GeV/c. 
Regardless of all the improvements in the comprehension of low-p^ production, the models are 
still unable to reproduce second order quantities such as final state particle correlations. We 
see that < > still increases as the multiplicity increases although not as fast. If we require 
Pt(Z) < 10 GeV/c, then HERWIG (without MPI) predicts that the < p T > decreases slightly 
as the multiplicity increases. This is because without MPI and without the high p? 'away-side' 
jet which is suppressed by requiring low pt{Z), large multiplicities come from events with a lot 
of initial-state radiation and the particles coming from initial-state radiation are 'soft', pythia 
tune AW describes the behavior of < pt > versus the multiplicity fairly well even when we 
select Pt{Z) < 10 GeV/c. This strongly suggests that MPI are playing an important role in 
both these processes. 
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7 Summary and Conclusions 

We are making good progress in understanding and modeling the softer physics. CDF tunes A 
and AW describe the data very well, although we still do not yet have a perfect fit to all the 
features of the CDF underlying event data. Future studies should focus on tuning the energy 
dependence for the event activity in the underlying event, which at the moment seems to be 
one of the least understood aspects of all the models. The underlying event is expected to be 
much more active in LHC and it is critical to have sensible underlying event models containing 
our best physical knowledge and intuition, tuned to all relevant available data. 
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Abstract 

We improve previously derived analytical estimates of hadronisation corrections to 
QCD jets at hadron colliders, firmly establishing at the two-loop level the link to the 
well-known power corrections to LEP event-shape variables. The results of this paper 
apply to jets defined in the kt and anti-kt algorithms but the general framework presented 
here holds also for other algorithms for which calculations are in progress. 

An understanding of QCD jets and their properties will be integral to the success of the 
LHC physics program. In particular one of the most important issues in current jet studies is 
the question of jet energy scale. The shift induced in jet energy by effects such as perturbative 
radiation and non-perturbative effects like hadronisation and the underlying event would con- 
tribute to a smearing of, for instance, mass peaks that may be a signal for new physics. Thus 
in order to choose optimal jet definitions that minimise such smearing one would need to know 
the dependence of the above effects on the experimental parameters like jet radius. Moreover, 
even in pure QCD studies such as the extraction of parton distribution functions (pdfs) and 
the strong coupling a s from jet observables like the inclusive jet cross-sections, a knowledge 
of the non-perturbative contribution is important to supplement perturbative calculations. A 
relatively small shift in the transverse momentum p t of a jet, induced by hadronisation, can 
result in a significant change in the inclusive jet spectrum since we are dealing with a quantity 
that has a steeply falling p t distribution. 

While it is traditional to study the hadronisation contribution via Monte Carlo models such 
as those in HERWIG and PYTHIA, it turns out that in cases like the jet energy there is additionally 
valuable analytical insight available [5 J . Analytical models based on renormalons pQ have in the 
past been met with great success in the description of LEP and HERA event-shape variables 
[7], but have not really been utilised outside that context. In Ref. [5] one such model (due to 
Dokshitzer and Webber [12]) was used to estimate hadronisation corrections to jet transverse 
momentum p t . The result found there was striking: hadronisation effects have a singular 1/R 
dependence on the jet radius R, at small R. This is in complete contrast to the contribution 
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from the underlying event which varies as R 2 . The knowledge of the R dependence of non- 
perturbative effects in conjunction with the \nR behaviour involved in perturbative estimates 
can then be used to arrive at conclusions about the optimal values of R to be used in diverse 
studies involving jets, as exemplified in Ref. [5]. 

While the computations of Ref. [5] indicate the dependence of hadronisation corrections 
and the underlying event on R, there remains the question of the overall magnitude of these 
effects. While for the underlying event one is reliant solely on Monte Carlo event generators 
to obtain the overall magnitude, for the hadronisation correction a tentative link was made in 
Ref. [5] between the magnitude of the 1/R correction and that of 1/Q corrections to LEP and 
HERA event shapes such as the thrust distribution (see [ZJ for a review). In order to definitely 
link the magnitude of jet hadronisation to that of event-shape power corrections one needs to 
carry out a calculation at the two-loop level rather than the simple one-loop estimate reported 
in [5]. The calculation for jets defined in the k t algorithm [131 E] is reported here while the 
corresponding result for jets in the anti-fct algorithm is already known [2J. Work on the other 
jet algorithms is currently in progress. 

1 The single-gluon result 

In the Dokshitzer- Webber model non-perturbative hadronisation corrections are associated to 
the emission of a soft gluon with transverse momentum k t ~ Aqcd- For the jet p t case we 
work out the change in transverse momentum 5p t induced by the emission of such a gluon and 
combined with the gluon emission probability (as given by perturbative QCD) this yields the 
average shift in p t induced by hadronisation: 



where k t , r] and <p respectively denote the transverse momentum, rapidity and azimuth of 
the emitted gluon with respect to the emitting hard parton (jet) and Cj is the colour factor 
associated to emission from the given hard parton. The reason we are able to single out a given 
hard parton initiating a high-p 4 jet in any hard process is essentially since the leading 1/R 
result stems from emission collinear to the triggered jet. It is thus possible to talk in terms of 
the single-jet limit ignoring the rest of the details of the hard process. 

The only non-perturbative ingredient that is involved above is the value of a s {k t ) at scales 
around or below Aqcd- If one makes the assumption of a universal infrared- finite coupling, 
which replaces the perturbative coupling that has an unphysical divergence at Aqcd, then 
one arrives at a prediction for the hadronisation correction (5p t ) h . Using the fact [5] that the 
Spt(k) is essentially the energy of the gluon emitted outside the jet we perform the integral 
over rapidity in Eq. (JTJ) to obtain the leading 1/R behaviour. The result is of the form cA/R, 
where c is a number obtained from the rapidity integral and A is the moment of the coupling 
&s{kt) over the infrared region (we refer the reader to Refs. [TQJ [TTJ for the precise details). 
Since the same coupling moment enters the predictions for event-shape variables we can take 
its value from data on event shapes and hence obtain a numerical prediction for the leading 
1/R hadronisation correction to jet p t . This was the method adopted in Ref. [5]. 




(1) 
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Figure 1: Gluon decay and one-loop corrections to single-gluon emission. 



Here we point out a limitation of the above approach [H] which is that while we have written 
down and used a running coupling a s (k t ), this quantity only emerges when one considers 
not just the emission of a single gluon but in fact gluon decay as well. To be precise an 
inclusive integration over gluon decay products is responsible for building up the quantity 
a s (k t ). Unfortunately, as is known for event-shape variables, our observable is sensitive to 
the precise details of gluon branching and hence one is not free to carry out such an inclusive 
integration. One must therefore return to the details of the gluon branching and identify the 
correction to the above inclusive approximation. The analysis at this level has already been 
carried out for event-shape variables [101 El El [8] and below we report on it for the jet p t case. 

2 Non-perturbative effects and gluon decay 

Now we consider the situation where the emitted gluon with k t ~ Aq CD is allowed to decay 
and at the same accuracy account for virtual corrections to single gluon emission, as depicted 
in Fig. CD 

At this two-loop level the change in p t (for a quark jet) can be expressed as: 



where a is a Sudakov variable, a s (0) is an ill-defined quantity which will cancel away subse- 
quently, x represents the virtual correction to gluon emission, dT 2 is the gluon decay phase-space 
and M 2 is the decay matrix element [101 [11] . We also denote by 5p t (ki, k 2 ) the change in p t 
due to correlated two-parton emission while Spt(k) is the corresponding single-gluon quantity. 
To correctly account for gluon branching one thus has to perform the above calculation, the 
details of which are reported in Ref. [I]. 

The analogous two-loop analysis for event-shape variables [101 [HI El E] revealed an initially 
surprising result - the two-loop correction simply provided a universal rescaling factor to the 
one-gluon result, which became known as the Milan factor. Its value for rif — 3 (which is 
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Figure 2: The role of clustering in determining whether a soft gluon ends up within or outside 
a hard jet. 



the number of flavours excited in the relevant soft region) was found to be M. = 1.49. Thus 
the ratio of corrections to two event shapes V\ and t>2 was merely the ratio of the one-loop 
coefficients computed previously [T2] : 



where 5v 1 ' denotes the non-perturbative single-gluon correction for v\ computed as discussed 
in the preceding section and likewise for V2- This remarkable result was understood to arise as 
a consequence of the fact that all the variables considered could be expressed as linear sums 
over the transverse momenta k ti of emissions, v = £\ kuCi, where the q are rapidity-dependent 
coefficients [TT] . 

In the case of jet p t this linear dependence is ruined by the non-trivial action of the jet 
algorithm in all cases except the case of jets defined in the anti-/c t algorithm |2J. The contribu- 
tion to the jet Spt of a given emission is found to be of the form kte 71 S out , where r\ denotes the 
rapidity with respect to the emitting hard jet and S out denotes the condition that the emission 
ends up outside the jet after the application of the jet algorithm. It should be immediately clear 
from this that in most current sensible jet algorithms (both of sequential recombination and 
cone type) the condition S out is non-trivial and introduces non-linearity in k t . For instance in 
the k t algorithm we can consider the situation in Fig. |2j where although one may have a soft 
parton separated by more than a certain distance R in rapidity and azimuth (denoted by the 
red gluon line) from a given hard parton, it may be clustered to another soft parton (denoted 
by the black gluon line) and hence swept into the final jet. This clustering depends on the 
k t of a soft parton relative to the other partons and hence the condition S out derived in [I] 
contains dependence on the k t of the soft partons, spoiling the simple linear dependence needed 
for universality. 

An exception to the above situation is to be found in the anti-fct algorithm for which the 
condition S out = Q{rf + (f> 2 — R 2 ) ensures that a given parton is outside the jet if its angular 
(r), 4>) separation from the hard jet is more than R. The linear dependence on k t is maintained 
and the Milan factor Ai = 1.49 is computed as for event shapes. 

For the k t algorithm we have carried out an equivalent calculation for the leading 1/R 
hadronisation correction (at small R) with the more complicated S out function involved there 
and we found the result M.^ t = 1.01 {rif = 3). Thus while at the level of the one-gluon studies 
of Ref. [5] the k t and anti-/c t algorithms received identical hadronisation corrections, a detailed 
analysis at the two-loop level breaks this equality. One finds that the ratio of hadronisation 
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corrections is then: 

= Hi -0.7. 
(fo^ ■ , 1-49 

Thus one expects somewhat smaller hadronisation corrections for the kt algorithm as com- 
pared to those for the anti-fct algorithm which is also borne out by the Monte Carlo studies with 
HERWIG and PYTHIA reported in Ref. [5]. We remind the reader that these conclusions apply 
only to the 1/R hadronisation corrections that would be dominant at small R and we neglect 
finite R corrections which need to be considered alongside the underlying event contribution 
which also has a regular R dependence ~ R 2 . 

3 Conclusions 

We have reported on a study of hadronisation corrections to jet pt or energy scale based on 
two-loop extensions of the one-gluon estimates reported in Ref. [5], with jets defined in the kt 
algorithm. Studies for other jet algorithms (SISCone [15] and Cambridge/Aachen [9]) are in 
progress. 
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